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The  presence  of  strain-induced  electric  fields  in 
lattice  mismatched  quantum  wells  (QWs)  grown  along  the  (111) 
growth  axis  gives  these  structures  advantages  over  those 
grown  along  the  (100)  growth  axis.  The  electric  field 
shifts  the  excitonic  resonance  to  lower  energy  (intrinsic 
Stark  effect)  and  decreases  the  optical  matrix  element.  This 
feature  is  absent  in  the  case  of  (100)  grown  layers  where 
the  piezoelectricity  along  this  axis  does  not  produce  any 
electric  field. 

In  this  work  we  present  the  theory  used  to  estimate  the 
strain-induced  electric  field  present  in  (111)  grown  In^Gai. 
3cAs/Alo.2Gao.8As  p-i-n  structures  for  x = 5%  and  10%.  The 
presence  of  this  field  was  proved  by  the  use  of  intensity- 
dependent  photoluminescence  measurements,  where  an  increase 
in  the  laser  light  intensity  used  in  the  PL  would 
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subsequently  increase  the  number  of  optically  created  e-h 
pairs  in  the  wells,  acting  to  screen  the  electric  field, 
causing  the  shift  of  the  excitonic  transition  energy  to  a 
higher  value. 

In  an  effort  to  further  investigate  the  effect  of  the 
strain,  a (111)  grown  sample  was  subjected  to  external 
mechanical  tensile  strain.  This  is  done  by  placing  the 
sample  on  a diaphragm  that  is  allowed  to  move  under  the 
application  of  a bias,  which  produced  a pulling  force  on  the 
sample,  i.e.  tensile  strain,  causing  a shift  in  the 
excitonic  peak  toward  a lower  transition  energy.  This  study 
thus  corroborated  the  existence  of  the  built-in  electric 
field  in  (111)  grown  samples. 

Finally,  the  nonlinear  refractive  index  for  the  samples 
grown  on  (lll)B  and  (100)  substrates  was  measured  using  a 
hybrid  Mach-Zender  interferometer  setup,  with  a TE  polarized 
laser  light  at  two  wavelengths  (980  and  992.2  nm) . An 
improvement  in  the  figure  of  merit  for  the  (111)  grown 
layers,  compared  to  those  grown  on  (100),  was  observed.  This 
finding  makes  the  use  of  (111)  grown  structures  as  nonlinear 
devices  more  advantageous  than  the  (100)  ones. 
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CHAPTER  1 
INTRODUCTION 


The  use  of  lattice  mismatched  material  has  been 
extensively  studied  for  structures  grown  on  (100)  GaAs 
substrates  for  their  use  in  optical  and  electrooptical 
devices  (mainly  having  InGaAs/GaAs  or  AlGaAs  MQW 
structure) [ 1-5] . More  recently  a new  effect  in  strained- 
layer  super  lattices  with  a [111]  growth  axis  has  been 
predicted.  It  is  the  presence  of  a large  internal  electric 
field  generated  by  the  piezoelectric  effect  [6].  Group  III-V 
semiconductors  are  piezoelectric;  thus,  strain  in  these 
materials  can  lead  to  electric  polarization  fields.  The 
magnitude  of  these  fields  can  be  very  large  exceeding  10^ 
V/cm  for  lattice  constant  mismatches  of  the  order  of  1%  [6]. 

The  strain  generated  electric  field  substantially 
changes  the  properties  of  (111)  strained  quantum  wells  and 
can  lead  to  large  nonlinear  optic  and  electrooptic  effects 
[7,8].  It  causes  the  tilt  of  the  energy  bands  and  leads  to  a 
red  shift  of  the  optical  transition  energies  [9,10]. 

In  this  dissertation  we  investigate  the  effects 
mentioned  above  for  compressively  strained  InGaAs /AlGaAs  MQW 
structures  grown  on  (lll)B  (As  terminated)  GaAs  substrates. 
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Chapter  2 explains  film  growth  procedure  on  (lll)B  GaAs 
substrates  using  a Molecular  Beam  Epitaxy  system.  The  growth 
of  GciAs,  AlGaAs  and  InGciAs  is  specially  emphasized.  We  also 
present  the  advantages  of  the  growth  along  this  axis  compared 
to  (100). 

Chapter  3 deals  with  the  theory  of  the  strain-induced 
electric  field  in  strained  MQWs  structures,  showing  the  main 
differences  between  (100)  and  (111)  grown  structures.  Then, 
two  structures  with  5%  and  10%  In  concentrations  are 
presented.  The  Kronig-Penney  model  is  introduced  in  order  to 
estimate  the  value  of  the  Cl-HHl  transition  energies  in  the 
samples,  keeping  in  mind  the  presence  of  the  electric  field 
in  the  (111)  grown  samples,  causing  the  transition  energy  in 
the  wells  to  shift  to  lower  energies.  This  effect  is 
interpreted  using  an  effective  well  width  model  proposed  by 
Miller  et  al.  [ 11 ] . 

Chapter  4 discusses  the  effect  of  the  increase  of 
optical  excitation  on  the  Cl-HHl  transition  energy.  The 
production  of  electron-hole  pairs  in  large  numbers  tends  to 
create  an  electric  field  which  screens  the  built-in  electric 
field,  and  in  turn  shifting  the  transition  energy  toward  a 
higher  energy  value  only  in  the  (111)  grown  structures.  This 
is  proven  by  the  use  of  a low  temperature  photoluminescence 
setup.  Furthermore,  the  effect  that  an  externally  applied 
strain  has  on  a (111)  grown  structure  is  measured  by 
monitoring  the  induced  energy  shift.  A simulation  is  run  to 
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determine  the  expected  change  in  the  layer  mismatch  and 
therefore  the  electric  field. 

Chapter  5 presents  the  means  used  to  measure  the  value 
of  the  nonlinear  refractive  index  in  all  the  structures.  It 
is  important  to  make  sure  that  the  compositions,  i.e. 
refractive  indices  of  the  layers  constituting  the  structures 
would  act  as  single  mode  slab  waveguides  at  the  wavelengths 
of  interest.  The  single  effective  oscillator  (SEO)  and  the 
multilayer  matrix  method  are  used  to  determine  the  refractive 
indices,  and  to  insure  a single  lateral  mode  operation, 
respectively.  Then,  a hybrid  Mach-Zender  setup  is  used  to 
obtain  the  values  of  U2  for  the  different  samples,  showing 
the  advantage  that  (111)  grown  samples  have  over  (100)  ones. 

Finally,  Chapter  6 summarizes  the  work  we  have 
undertaken  in  the  course  of  this  research,  giving  some 
predictions  and  recommendations  for  further  studies. 


CHAPTER  2 

FILM  GROWTH  ON  (lll)B  GaAS 
2.1  Introduction 


III-V  compound  semiconductors  have  been  grown 
traditionally  on  (100)  substrates  because  of  the  wide  range 
of  growth  conditions  which  result  in  good  epitaxial  layer 
quality  on  the  this  crystalline  orientation.  However, 
fundamental  material  properties,  growth  mechanisms,  surface 
kinetics,  and  quantum  well  properties  can  be  changed  or 
improved  as  compared  to  those  on  (100)[9,12,13].  It  has  been 
theoretically  predicted  and  observed  that  coherently  strained 
(110)  and  (111)  layers  will  possess  a large  built-in  in-plane 
or  perpendicular-to-plane  piezoelectric  field,  which  leads  to 
a new  class  of  self  electro-optic  devices.  Quantum  wells 
grown  on  (111)  possess  an  enhanced  optical  transition  [14], 
and  a low-threshold  (111)  QW  laser  has  been  made  [12]. 

2.2  Growth  on  (111)  Oriented  Substrates 
2.2.1  Introduction 


It  is  rather  difficult  to  grow  high  quality  materials 
on  (111)  substrates  [15].  Either  forced  two-dimensional 
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layer-by-layer  migration  enhanced  epitaxy  or  a very  high 
growth  temperature  (700°C)  are  required  to  suppress 
microtwins  and  facet  growth  on  (111).  However,  Chin  et 
al.  [16-18]  were  able  to  grow  high  quality  layers  of  GaAs, 
AlGaAs  and  InGciAs  on  (lll)B  GaAs. 

To  achieve  this  goal,  Chin  et  al. [16-18]  have  examined  a 
variety  of  miscuts  on  (lll)B  GaAs.  Few  degrees  of  miscut  can 
provide  growth  ledges  to  achieve  two-dimensional  step  flow 
epitaxy.  It  was  found  that  if  (lll)B  was  miscut  toward  <211> 
instead  of  <100>  and  <011>,  a wide  range  of  angles  (l°-3°)  can 
generate  a smooth  growth.  Although  it  was  also  found  that 
0.5°  miscut  toward  <10 0>  can  achieve  a very  high  quality 
material  [12],  it  was  not  used  in  the  growth  of  the  proposed 
structures  because  an  exact  0.5°  miscut  is  extremely 
difficult  to  control.  A proper  choice  of  miscut  is  the  key 
in  achieving  a high  quality  material  on  (111),  over  a wide 
range  of  growth  conditions.  To  discover  which  of  the  miscut 
angles  is  best  for  growth,  GaAs  and  AlGaAs  were  grown  on 
(lll)B  GaAs. 

2.2.2  Molecular  Beam  Epitaxy  (MBE)  Growth 

The  essential  elements  of  a MBE  system  [19]  are  shown  in 
Fig.  2.1.  It  is  apparent  from  this  illustration  that  each 
MBE  arrangement  may  be  divided  into  three  zones  where  three 
different  physical  phenomena  take  place  [20].  The  first  zone 
is  the  generation  zone  of  the  molecular  beams . Next  is  the 
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zone  where  the  beams  from  different  sources  intersect  each 
other  and  the  vaporized  elements  mix  together  creating  a very 
special  gas  phase  contacting  the  substrate  area.  This  area, 
where  the  crystallization  processes  take  place,  can  be 
regarded  as  the  third  zone  of  the  MBE  physical  system. 

The  molecular  beams  are  generated  in  the  first  zone 
under  UHV  conditions  from  sources  of  the  Knudsen-ef fusion- 
cell type  [20],  whose  temperatures  are  accurately  controlled. 
Conventional  temperature  control,  based  on  high  performance 
proportional-integral-derivative  (PID)  controllers  and 
thermocouple  feedback,  enables  a flux  stability  of  better 
than  ±1%  [21].  By  choosing  appropriate  cell  and  substrate 
temperatures,  epitaxial  films  of  the  desired  chemical 
composition  can  be  obtained.  Accurately  selected  and 
subsequently  precisely  controlled  temperatures  of  the 
substrates  and  each  of  the  sources  of  the  constituent  beams 
have  thus  a different  effect  upon  the  growth  process. 

The  second  zone  in  the  MBE  vacuum  reactor  is  the  mixing 
zone,  where  the  molecular  beams  intersect  each  other.  Little 
is  known  at  present  about  the  physical  phenomena  occurring  in 
this  zone.  This  results  from  the  fact  that  usually  the  mean 
free  path  of  the  molecules  belonging  to  the  intersecting 
beams  is  so  long  that  no  collisions  and  no  other  interactions 
between  the  molecules  of  different  species  occur  there. 

Epitaxial  growth  in  MBE  is  realized  in  the  third  zone, 
i.e.  on  the  substrate  surface.  A series  of  surface  processes 
is  involved  in  MBE  growth;  however,  the  following  are  the 
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Fig  2.1  Schematic  illustration  of  the  essential 
parts  of  an  MEE  growth  system  [19]. 
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most  i mportant  [22]  : 

i)  Adsorption  of  the  constituent  atoms  or  molecules 
impinging  on  the  substrate  surface. 

ii)  Surface  migration  and  dissociation  of  the  adsorbed 
molecules . 

iii)  Incorporation  of  the  constituent  atoms  into  the  crystal 
lattice  of  the  substrate  or  the  epilayer  already  grown. 

iv)  Thermal  desorption  of  the  species  not  incorporated  into 
the  crystal  lattice. 

These  processes  are  schematically  illustrated  in  Fig  2.2 
[23].  The  substrate  crystal  surface  is  divided  into  so- 
called  crystal  sites  on  which  the  impinging  molecules  or 
atoms  may  interact.  Each  crystal  part  is  a small  part  of  the 
crystal  surface  characterized  by  its  individual  chemical 
activity.  A site  can  be  created  by  a dangling  bond,  vacancy, 
step  edge,  etc.  [24]. 

2.2.3  Growth  of  GaAs,  AlGaAs,  and  InGaAs 

A Varian  Modular  GEN  I I Molecular  Beam  Epitaxial  growth 
system  was  used  to  grow  GaAs  and  AlGaAs  on  (lll)B  and  (100) 
substrates  simultaneously.  This  system  was  carefully  baked 
and  outgassed  to  reduce  background  impurities.  After  a few 
days  of  bakeout  and  cooling,  the  chamber  pressure  is  lower 
than  1x10“^^  Torr.  The  careful  outgas  procedure  and  low 
ambient  pressures  of  CO,  H2O,  and  O2  molecules  are  important 
to  achieve  high  quality  growth  with  low  concentrations  of 
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Fig  2.2  Schematic  illustration  of  the  surface  processes 
occurring  during  film  growth  by  MBE  [22]. 
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nonradiative  recombination  centers. 

The  quality  of  the  grown  material  was  characterized  by 
the  use  of  low  temperature  photoluminescence  (PL) 
measurements.  Chin  et  al.[12]  found  that  the  principal 
spectral  features  of  the  GaAs  layer  grown  on  (100)  are 
neutral  donor  bound  (D°,X)  and  acceptor  bound  (A°,X)  exciton 
recombination  lines  and  free  electron-acceptor  (e,A°)  and 
donor  acceptor  (Do,A°)  transitions.  However,  only  the  neutral 
donor  bound  exciton  transition  was  observed  for  the  GaAs 
layer  grown  on  (lll)B. 

The  acceptor  related  transition  for  (100)  is  due  to 
carbon  incorporated  into  the  GaAs  layer  from  MBE  background 
carbon-related  species  such  as  CO,  CH4,  etc.  The  donor 
related  transition  is  due  to  sulfur  in  high  quality  MBE 
grown  GaAs. 

Low  carbon  acceptor  incorporation  in  (lll)B  GaAs  can  be 
understood  based  on  the  number  of  dangling  bonds.  While  the 
(100)  surface  always  has  two  double  dangling  bonds  available 
for  both  Ga  and  As  surface  adatoms,  the  (lll)B  surface  has 
either  one  or  three  dangling  bonds  available  for  Ga  and  As, 
respectively  (see  Fig.  2.3).  The  large  number  of  dangling 
bonds  available  to  As  on  (lll)B  increase  the  As  sticking 
coefficient  and  reduce  surface  lifetime  for  As  adatoms.  Thus 
the  probability  for  carbon  incorporation  is  much  reduced. 
This  is  one  of  the  advantages  of  material  growth  on  (lll)B 
since  low  carbon  incorporation  in  GaAs  improves  the  material 
transport  properties. 
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Next,  Chin  et  al.[18],  compared  the  growth  of  2iim 
Alo.22Gao.78As  on  (100)  and  (lll)B  GciAs.  Growth  was  carried 
out  simultaneously  on  1°,  2°  and  3°  misoriented  (lll)B  and 
0°(100)  GaAs  substrates  at  630  and  670°C,  respectively. 

The  PL  spectra  for  the  growth  temperature  of  630°C  are 
shown  in  Fig.  2.4.  A full  width  at  half-maximum  (FWHM)  of 
4.4  meV  for  (100)  AlGaAs  indicates  excellent  quality.  The 
high  energy  peak  at  1.837  eV  is  due  to  bound  excitons.  The 
broad  peak  centered  at  1.818  eV  is  a combination  of  electron- 
to-neutral-carbon-acceptor  and  neutral-donor-to-neutral- 
acceptor  recombination. 

The  (lll)B  AlGaAs  layers  show  broad  peaks  without  any 
resolvable  fine  structures.  The  peak  energies  decrease 
monotonically  from  1.837  eV  for  0°(100)  to  1.744  eV  for  the 
3°(111)B  layer  (a  shift  of  93  meV) . The  1°  and  3°  misoriented 
layers  show  single  peaks  almost  10  times  as  intense  as  that 
of  the  (100),  whereas  the  2°  misoriented  layer  shows  two 
peaks  half  as  intense  as  those  of  1°  and  3°  misorientations . 
The  fact  that  the  PL  spectra  are  not  due  to  impurity  related 
transitions  is  evidenced  by  low  impurity  content  of  the 
layers  indicated  by  Hall  mobility  measurements.  The  FWHM  for 
1°  and  3°  misorientations,  however,  increases  to  17.5  and 
20.2  meV,  respectively.  The  half-widths  at  half-maximum 
(HWHM)  for  2°  misorientation  peaks  are  9.5  and  20.2  meV. 

Figure  2.5  shows  the  PL  spectra  of  the  AlGaAs  grown  at 
670OC.  The  energy  of  the  peak  of  the  0°(100)  layer  is  about  8 
meV  higher  than  that  grown  at  630°C.  This  shift  in  energy  is 
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due  to  higher  Ga  desorption  (re-evaporation  of  surface  Ga 
adatoms)  at  this  growth  temperature,  with  a slight  increase 
in  the  linewidth  to  5.0  meV.  The  peak  intensities  for  (lll)B 
show  little  change  from  that  grown  on  0°(100)  substrate 
although  a decrease  in  the  linewidth  is  observed  from  a 16.4 
meV  for  1°  to  8.8  meV  for  3°  misorientation.  All  the  peak 
energies  lie  within  9 meV  of  that  for  the  (100)  layer  and  the 
maximum  increase  in  the  integrated  intensity  is  by  a factor 
of  2.7  only.  The  absence  of  any  significant  orientation 
dependence  of  the  PL  for  the  layers  grown  at  67  0°C  is 
believed  to  be  due  to  high  surface  migration  of  both  A1  and 
Ga  adatoms. 

The  relative  PL  integrated  intensity  and  peak  energy 
shifts  are  summarized  in  Figs.  2.6a  and  2.6b  , respectively. 
An  order  of  magnitude  intensity  enhancement  in  (100)  AlGaAs 
grown  at  670°C  compared  to  630°C  is  due  to  reduction  in  the 
defect  related  nonradiative  recombination  centers.  This  is 
expected  from  higher  surface  migration  velocities  of  adatoms 
and  lower  oxygen  incorporation  at  higher  growth  temperature. 

On  the  contrary,  a decrease  in  the  PL  integrated 
intensities  for  (111)  AlGaAs  as  the  growth  temperature  is 
increased  cannot  be  explained  by  the  above  argument. 

Furthermore,  to  characterize  InGaAs  growth  on  (lll)B 
GaAs,  Chin  et  al.  [17]  have  grown  a 20  A Ino.21Gao.79As  QW  at 
550°  C,  for  which  the  linewidth  of  the  PL  spectrum  at  2.2  K is 
extremely  small  (only  3.0  meV),  showing  the  high  quality  of 
the  growth. 
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Fig  2.3  A (111)  surface  of  a diamond  type  crystal  [24]. 
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Fig  2.4  13  K PL  of  2.0  ^thick  AIq. 22 Gao. 79 As  grown  at  630°C 
[18]. 
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Fig.  2.5  13  K PL  of  2.0  /mi-thick  AIq. 22 Gag. 79 As  grown  at  670°C 
[18]. 
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The  InGaAs/AlGaAs  structures  we  investigated  were  grown 
on  2°  misoriented  undoped  (lll)B  and  (100)  GaAs  substrates. 
The  AlGaAs  and  the  InGciAs  layers  were  grown  at  substrate 
temperatures  of  670  and  550  *C  respectively.  The  reason  for 
choosing  2°  misoriented  (lll)B  GaAs  substrate  and  a growth 
temperature  of  670 *C  can  be  clearly  seen  in  Figs  2.4  and  2.5. 
The  PL  peak  intensity  of  the  AlGaAs  layer  grown  at  670 “C  and 
on  2°  misorientation  (lll)B  GaAs  is  much  higher  than  all 
other  misorientations  at  both  growth  temperatures,  and  its 
FWHM  is  also  smaller  than  that  of  the  3°  misoriented 


substrate . 
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Fig.  2.6  Substrate  orientation  dependence  of  [18] 

(a)  EL  integrated  intensity. 

(b)  PL  peak  energy 

for  Alo.22Gao.79As  grown  at  630  and  670°C . 


CHAPTER  3 

CHARACTERISTICS  OF  STRAINED  MQW  STRUCTURES 
GROWN  ON  (lll)B  GaAs  SUBSTRATES 


3.1  Introduction 


Strained  MQW  structures  grown  on  (lll)B  GaAs  substrates 
differ  significantly  from  their  (100)  counterparts.  Layers 
grown  along  the  (111)  direction  exhibit  an  increase  in  the 
value  of  the  heavy  hole  effective  mass  conversely  rendering  a 
lower  value  for  the  HHl  transition  in  the  wells  and 
decreasing  the  Cl-HHl  transition  energy. 

Moreover,  the  presence  of  strain  due  to  the  growth  of 
lattice  mismatched  layers (namely  InGaAs  on  (lll)B  GaAs) 
produce  a strain-induced  electric  field  in  these  layers  due 
the  piezoelectric  nature  of  the  (111)  grown  layers,  yielding 
a shift  toward  a lower  Cl-HHl  transition  energy  due  the 
Quantum  Confined  Stark  Effect  (QCSE). 

On  the  other  hand,  due  to  the  hydrostatic  and  shear 
strain  present  in  the  lattice  mismatched  layers,  the 
conduction  to  the  valence  band  transition  energy  will  go  up 
significantly  in  the  case  of  (111)  grown  layers  compared  to 
that  of  (100)  grown  layers. 

After  discussing  the  above  mentioned  characteristics,  we 
introduce  our  suggested  p-i-n  structures  grown  on  (lll)B  and 
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(100)  GaAs.  Then,  the  characteristics  of  the  layers  grown 
along  these  orientations  will  be  discussed  showing  the  effect 
that  this  growth  orientation  has  on  the  performance  of  these 
layers.  Finally,  we  explain  the  reasons  behind  choosing 
these  structures,  showing  the  advantage  that  our  structures 
has  over  other  structiures  introduced  by  previous  researchers, 
and  their  possible  use  as  optical  devices,  mainly  Self- 
Electrooptic  Effect  Devices  (SEEDs). 

3.2  Characteristics  of  strained  Structures  on  Qll^B 

Substrates 


3.2.1  Effectivp  Mass  Changes  due  to  the  Growth  Orientation 


The  effective  mass  change  affects  the  quantization  shift 
in  the  grown  layers  and  is  also  understood  to  influence  the 
optical  matrix  element  and  the  density  of  states.  It  was 
shown  [25-27  that  the  heavy  hole  effective  mass  along  the 
growth  axis  changes  considerably  between  the  (100)  and  (111) 
growth  orientations.  The  heavy  hole  and  light  hole 
effective  masses  along  the  growth  axis  for  the  (100)  and 
(111)  oriented  structures  are  expressed  as  [14] 


m 
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Ih_ 
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Yi±2y2 
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Wl 
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(3.1) 
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where  Yl»Y2  2ndY3  are  the  Luttinger-Kohn  [28,29]  parameters 
and  mo  is  the  free  electron  mass. 

Using  the  values  given  in  Table  3.1,  for  GaAs, 
for  exanple,  is  0.917mo  vdiile  mhh(lOO)  for  the  same  material 
is  0.353mof  where  mo  is  the  free  electron  mass.  Since  the 
(111)  oriented  heavy  hole  mass  is  2.6  times  higher  than  that 
of  the  (100)  orientation,  the  quantized  energies  in  the  QWs 
for  heavy  hole  subbands  in  the  (111)  structure  therefore  are 
smaller.  As  a result,  even  higher  order  transitions  between 
electron  and  heavy  hole  subbands  are  expected  to  be  observed 
in  transmission  measurements  for  the  (111)  oriented  QWs.  On 
the  other  hand,  the  electron  effective  mass  is  the  same  for 
both  (111)  and  (100)  grown  samples,  while  the  light  hole 
effective  masses  vary  only  by  a small  amount  (mih(lll)=0.07mo 
and  iHih(100)=0»08mo  for  GaAs).  Thus,  the  light  hole  related 
transitions  energies  are  considered  to  be  almost  the  same  in 
(100)  and  (111)  oriented  QWs. 


Table  3.1 

Luttinger-Kohn  Parameters 


Y, 

GaAs 

7.65 

2.41 

3.28 

InAs 

19.67 

8.37 

9.29 

AlAS 

4.04 

0.78 

1.57 

3.2.2  Piezoelectric  Effec± 
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Unlike  (100)  oriented  strained  layers,  (111)  strained 
layers  of  piezoelectric  active  materials,  such  that  zinc- 
blende  confounds,  have  built-in  longitudinal  electric  fields 
[ 6 ] . Piezoelectric  active  crystals  are  a class  of  crystals 
which  produce  a macroscopic  electric  polarization  when 
stressed.  Because  electric  polarization  is  a vector, 
crystals  \diich  have  a center  of  inversion  symmetry  cannot  be 
polarized  by  stress  and  therefore  cannot  be  piezoelectric 
[30].  Zinc-blende  crystals  lack  inversion  symmetry  and  are 
known  to  be  piezoelectrically  active  [31].  Stressing  a 
crystal  produces  a microscopic  displacement,  or  strain,  of 
the  atoms  from  their  equilibrium  positions.  If  the  strain 
results  in  an  unequal  motion  of  the  positive  charges  with 
respect  to  the  negative  charges,  then  an  electric 
polarization  will  be  produced.  Misfit  strain,  generated  in  a 
crystal  during  pseudomorphic  epitaxy,  is  accommodated  by  the 
competing  process  of  elastic  deformation  and  misfit 
dislocation  formation,  both  of  which  act  to  minimize  the 
strain  energy  of  the  system  [32].  In  sufficiently  thin 
overlayers,  vdiere  misfit  dislocations  are  energetically  not 
favored,  the  overlayer  is  constrained  to  grow  in  registry 
with  the  negligibly  strained  substrate.  In  this  way,  the 
layer  is  subjected  to  a biaxial  stress,  in  the  plane  of  the 
growth,  which  induces  an  elastic  deformation,  leaving  the 
overlayer  tetragonally  distorted  [33].  Although  both  (100) 
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and  (111)  oriented  strained  crystals  cire  grown  under  biaxial 
stress  and  experience  a longitudinal  tetragoneil  distortion, 
only  the  (111)  growth  axis  strained  layers  produce  a 
macroscopic  polarization.  This  polarization  is  given  by 

Pi  = ©14  6jk  (3.3) 

where  is  the  induced  polarization,  03^4  is  the 
piezoelectric  constant  [34],  and  ejj^  is  the  symmetrized 

strain  conponent  [35],  However,  diagonal  conponents  of  the 
strain  matrix  (e.g.  do  not  induce  a polarization  (i.e., 
ej^j^  = 0)  in  the  zinc-blende  materials  [31].  The  sign  of  the 
polcirization  vector  in  a material  depends  on  whether  it  has 
the  larger  or  smaller  lattice  constant  and  on  the  sign  of  the 
piezoelectric  coefficient.  The  common  III-V  semiconductors 
have  a negative  piezoelectric  coefficient  [35].  A III-V 
semiconductor  with  the  larger  lattice  constant  in  a strained 
layer  superlattice  will  have  the  polarization  vector  pointing 
from  the  A (cation)  to  the  B (anion)  face.  A III-V 
semiconductor  with  the  smaller  lattice  constant  will  have  the 
polarization  vector  pointing  frcan  the  B to  the  A face.  Thus, 
for  a strained-layer  superlattice  with  [111] -growth  axis  of 
two  III-V  semiconductors,  polarization  vectors  along  the 
growth  axis,  and  of  opposite  sign  in  the  two  materials ( InGaAs 
and  AlGaAs  for  the  case  studied  here ) , are  induced  by  the 
piezoelectric  effect.  Furthermore,  for  the  case  we  are  about 
to  present  where  the  InGaAs  layer  is  under  biaxial 
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con^)ression,  the  strain-induced  electric  field  will  point 
toward  the  substrate  for  a (lll)B  (i.e.  As  terminated) 
substrate  (see  Fig.  3.1). 

Because  the  three  off-diagonal  stress  con^nents  are 
equal,  the  components  of  the  polarization  vector  are  equal 
and  the  polarization  vector  is  in  the  [111] -growth 
direction.  Thus,  there  are  four  unknowns:  a diagonal  cind  an 
off-diagonal  strain  component  in  each  material.  Requiring 
that  the  strain-distorted  lattice-translation  vectors  of  the 
two  materials  have  equal  projections  on  the  plane  normal  to 
the  growth  direction  gives 

d + e«)-^a*  = (l  + E^)-^a‘’  (3.4) 


where  a^  is  the  lattice  constant  of  material  i.  The 
mechanical  energy  density  of  the  superlattice  is  minimized 
with  respect  to  the  strain  components  subject  to  the 
condition  of  Eq.  (3.4).  This  procedure  gives 
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(3.5d) 


(3.5e) 


(3.5f) 


vrtiere  hi  is  the  layer  thickness  of  matericd  i and  the  C's  are 
the  elastic  coefficients  of  the  materials  a and  b.  For  (111) 
grown  layers,  the  diagonal  stress  ccnponents  ( s*xx/  ^yy,  eL)  do 

not  induce  any  internal  electric  field,  only  off-diagonal 
stress  conponents  ( ejjy,  6^,  E^)  need  to  be  considered. 

The  strain  induced  electric  polarization  Pi  will  lead  to 
electric  field  Ei  given  by 


Di  - £q  X Ei  + Pi 


(3.6) 


where  x is  the  susceptibility  and  Di  is  the  electric 
displacenent.  If  there  are  no  external  charges,  Di  vanishes, 
and  the  electric  field  reduces  to 


Ei  = 


^14gjk 

Eo£ 


(3.7) 


vdiere  E = 1 + x is  the  low  frequency  dielectric  constant. 

The  strain  induced-field  tends  to  shift  the  electron 
and  hole  quantized  energies  by  QCSE,  thus  reducing  the 
effective  band  gap  [14]. 
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Strain-induced  Electric  Field 


3.2.3  Band  Structure  Modification  due  to  st.rain  Effect 
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For  conditions  that  the  thickness  of  the  layer  is  less 
than  the  critical  thickness  [36-38],  the  streiin  is  confined 
in  the  layer  without  the  formation  of  dislocations 
(pseudomorphic  growth).  In  this  case,  the  lattice  constant 
of  the  epitaxial  layers  in  the  plane  perpendicular  to  the 
growth  axis  is  equal  to  that  of  the  substrate.  As  a first- 
order  perturbation  of  the  strain  Hamiltonian,  the  strain 
induced  band  structure  modification  in  the  (111)  grown  layers 
can  be  written  as 

AEc  = 3 3c  £\x  (3.8) 

for  the  conduction  band  and 

AEhh  = 3 3v  £xx  ± d E,y  (3.9) 

Ih 

for  the  heavy  hole  band  (+)  and  the  light  hole  band  (-), 
respectively  [14],  a^  and  a^  are  the  conduction  and  the 

valence  bands  hydrostatic  deformation  potentials,  and  d is 
the  shear  deformation  potential.  In  Eqs.  3.8  and  3.9,  energy 
positive  directions  are  set  to  be  opposite  between  conduction 
and  valence  band.  It  was  shown  [33]  that,  for  ccxrpressively 
strained  layers  (0.5  % strain)  grown  on  (100)  and  (111) 
oriented  InP  substrates,  the  valence  band  splitting  due  to 
shear  strain  was  very  similar,  whereas  the  bandgap  change 
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due  to  hydrostatic  strain  was  much  larger  for  (111)  than  for 
(100).  The  Cl-HHl  transition  for  the  (111)  grown  layers  was 
33  meV  higher  than  the  unstrained  case  \diile  for  (100)  it  was 
only  17  meV  higher  (see  Fig.  3.2). 

3.3  A Study  of  rilDB  Grown  Strained  MCW  Structures 
3.3.1  Introduction 

The  structures  investigated  are  p-i-n  diodes,  with  the 
wells  imbedded  in  the  intrinsic  region  of  the  diode,  having 
different  In  conpositions  and  therefore  different  strcdn  and 
built-in  strain-induced  electric  fields. 

The  p-i-n  structures  were  grown  simultaneously  by  MHE  on 
(lll)B  and  (100)  GaAs  substrates.  A schematic  of  the 
structures  is  shown  in  Fig. 3. 3 (a)  and  (b).  The  first 
structure  (SI)  is  composed  of  a 0.8  pm  n"*"  GaAs  layer, 
followed  by  a 2.5  fun  n"*"  AIq  2 Gag^gAs  layer,  and  undoped  250 
nm  AIq  2^]^  Gag  g^As  layer,  then  the  MQW  region  (10  periods,  15 
nm  Alg  2 Gag  gAs  / 10  nm  Ing  j^Gag  gAs),  then  an  undoped  250  nm 
Alg  11  Gag  ggAs  layer,  and  a 1 fun  p‘*‘  Alg  2Gag  gAs  layer  capped 

with  a 50  nm  p"^  GaAs.  The  structures  were  designed  to  work 
as  single  mode  slab  waveguides  for  wavelengths  0.8-1. 3 pm. 

The  second  structure  (S2)  is  similar  to  the  first  but  has 
Alg  ivGag  ggAs  instead  of  Alg  Gag  ggAs  cladding  layer  and 

also  iHg  gjGag  95AS  Instoad  of  Ing  ^Gag  gAs  in  the  well. 
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InGaAs  on  AIGaAs 


Strained  on  ( 111 ) Strained  on  ( 1 00) 


Fig  3.2  Different  energy  band  shifts  for  (100)  and  (111). 
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Fig.  3.3  Structure  of  proposed  p-i-n  devices. 


3.3.2  Kroniq-PenneY  Model 
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A theoretical  model  based  on  the  Kronig-Penney  method 
[39]  was  used  to  ccd.culate  the  Cl-HHl  MQW  treinsitions  in  the 
above  mentioned  structures. 

The  heavy  hole  effective  masses  of  the  GaAs,  AlAs  and 
InAs  are  calculated  using  equation  3.2.  Then,  the  mhh  of 
InGaAs  and  AlGaAs  are  obtained  using  the  following 
interpolation 

mhhAiGaAs  “ X ®hhAlAs  (1“X)  mhhCaAs  (3.10) 
n^ihinGaAs  = Y rohhinAa  + (l~y)  n^hhcaAs  (3.11) 


similarly,  the  electron  effective  masses  is  of  AlGaAs  and 
InGaAs  are  given  by 

mcAlGaAs  ~ nt*(0.0665  + 0.083  x)  (3.12) 

incinGaAs  = nb*(y  0.023  + (1-y)  0.067  (3.13) 

vdiere  x and  y are  the  Al  and  In  fractions,  respectively  and 
itio  is  the  free  electron  mass. 

In  the  same  manner,  the  lattice  constants  for  AlGaAs  and 
InGaAs  are  obtained  by  interpolating  the  values  of  the 
lattice  constants  of  GaAs,  AlAs  and  InAs  mentioned  in  Table 

3.2. 
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The  bandgap  energies  of  the  AlGaAs  barriers  and  the 
InGaAs  wells  are  determined  using  the  following  equations  : 

S^AlGaAs  = 1.424  + 1.266  x + 0.26  (3.14) 

EgmcsaAs  = 0.36  + 0.505  (1-y)  + 0.555(l-y)2  (3.15) 


Eg  - EgAlGaAs  - EginGaAs  (3.16) 

The  relation  between  the  conduction  band  discontinuity 
AEc  and  Eg  is  for  InGaAs /Alo.2Gao.8As  quantum  well  structures 
given  by  [40] 

AEc  = 0.6  Eg  (3.17) 


vrtiile  the  valence  band  discontinuity  is  given  by 


AEv  = Eg  - AEc  (3.18) 

The  quantization  energy,  Ed,  of  the  first  electron  level 
in  the  conduction  band  is  equal  to  the  energy  for  vrtiich 


sinh  yL.  sinaL,  + cosh  yLuCosaL  =-l 

D Z ' D Z 


(3.19) 
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Y = 


(AEc-Eci)mCjnGaAs 


(3.21) 


v^ere  Lb  and  Lz  are  the  barrier  and  well  width,  respectively. 

Similarly,  the  value  of  Evl,  the  first  electron  level  in 
the  valence  band  is  acquired  using  the  equations  above,  after 
exchanging  mcmcaAs  by  mhhmGaAs/  Eel  by  Evl  and  AEc  by  AEv 
Fig.  3.4  shows  the  Cl-HHl  transition  without  the  presence  of 
streiin  or  electric  field. 


Table  3.2 

Parameters  of  GaAs,  AlAs,  inAs  and  InGaAs 


ao 

(^) 

Cll*10^ 
dynes /cm2 

Ci2*10^^ 
dynes /cm2 

044*10^1 
dynes /cm2 

ei4 

GaAs 

5.6532 

11.810 

5.32 

5.940 

-0.1600 

AlAS 

5.6611 

12.020 

5.70 

5.870 

-0.2250 

InAs 

6.0584 

8.330 

4.53 

3.960 

-0.0450 

InGaAs 
(5%  In) 

5.6734 

11.636 

5.28 

5.841 

-0.15425 

InGaAs 
(10%  In) 

5.69372 

11.462 

5.241 

5.742 

-0.1485 
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Fig.  3.4  Energy  band  structure  of  a QW  with 
no  electric  field  present. 


3.3.3  s-fcrain  Effect  on  Cl-HHl  Transition  Energy 
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would  move  the  electron  to  heavy  hole  transition  to  a higher 
energy.  For  the  (lll)B  grown  layers,  this  shift  can  be 
determined  by  replacing  the  values  of  ac,  av  and  d for  each 
InGaAs  composition  in  equations  3.8-9.  For  the  same 
structures  grown  on 

(100)  GaAs  substrates,  the  hydrostatic  and  shear  strains  will 
cause  a similar  shift  but  by  a different  amount.  The  change 
in  this  case  is  given  by  [41] 


AE 


Hi 


+ b 


+ 2 C 


'u 


(3.22) 


Where  is  the  shift  of  the  heavy  hole  valence  band,  £ is 

the  biaxial  strain  parallel  to  [100]  and  [010]  directions, 
taken  here  to  be  the  mismatch  Aa  given  by 


Aa  - ( ainGaAs“^AlGaAs ) /SAlGaAs 


(3.23) 


a and  b are  the  hydrostatic  and  the  shear  deformation 
potentials,  respectively.  The  values  of  a and  b are  given  in 
Table  3.3  for  GaAs  and  InAs.  These  values  for  InGaAs  are 
calculated  by  interpolating  their  respective  binary  values. 

Conparing  the  amount  of  shift  between  the  (100)  and  the 
(lll)B  grown  structures  will  show  a higher  value  in  the  case 
of  the  latter.  This  presents  an  advantage  in  using  the 
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(lll)B  grown  structures  yielding  higher  change  in  the 
transition  energies  for  the  same  Aa. 


Table  3.3 

a,  b,  ac,  av  and  d Parameters  of  GaAs  and  InAs 


GaAs 

InAs 

a(eV) 

2.7 

2.5 

b(eV) 

-1.7 

-1.8 

ao(eV) 

-18.3 

-11.7 

av(eV) 

-4.8 

-5.2 

d(eV) 

-4.4 

-4.6 

3.3.4  Electric  Field  Effect  on  CWs 

The  presence  of  strain-induced  electric  field  in  the 
wells  will  cause  a red  shift  in  the  transition  energies. 
This  can  be  explained  by  the  Quantum-Confined  Stark  Effect 
(QCSE).  This  effect  can  be  explained  as  follows.  An 
electric  field  perpendicular  to  the  quantum  well  layers  pulls 
the  electrons  and  holes  towards  opposite  sides  of  the  layers 
resulting  in  an  overall  net  reduction  in  energy  of  an 
electron-hole  pair  and  a corresponding  Stark  Shift  in  the 
exciton  absorption  (see  Fig.  3.5).  The  walls  of  the  quantum 
well  inpede  the  electron  and  hole  from  tunneling  out  of  the 
>^11  in  rapid  field  ionization.  Because  the  well  is  narrow 
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(in  our  case  lOOA)  compared  to  the  three  dimensional  (3D) 
exciton  size  (e.g.  300A),  the  electron-hole  interaction, 
although  weakened  by  the  separation  of  electron  and  hole,  is 
still  strong  and  well  defined  excitonic  states  can  still 
exist.  Kajikawa  et  al.[42]  have  showed  an  increase  in  the 
rate  of  energy  shift  of  the  Cl-HHl  exciton  under  the 
application  of  an  electric  field  for  (111)  C3Ws  over  that  for 
(100)  QWS  having  the  same  width. 

To  determine  to  vdiat  extent  the  presence  of  the  strain- 

induced  electric  field  would  affect  the  Cl-HHl  transition,  we 

will  use  the  effective  infinite  well  model  reported  by  Miller 

et  al . [11].  This  model  uses  the  infinite  well  model  with 

effective  well  widths  assigned  separately  for  the  electron 

and  the  heavy  hole,  instead  of  the  actual  well  size,  in  order 

to  obtain  the  correct  zero-field  energies.  The  exact 

solution  for  a particle  in  an  infinite  well  in  the  presence 

of  a uniform  static  field  is  explained  as  follows.  It  was 

found  [43-45]  that  if  the  particle  is  restrained  to  stay 

within  an  infinitely  deep  well,  an  exact  solution  can  be 

obtained  for  any  value  of  the  static  field.  This  is  the 

consequence  of  the  boundary  conditions  inposed  upon  the  wave 
function  ^(z),  i.e.,  ^(±i  Lz)  = 0 for  any  field. 

The  Schrodinger  equation  can  be  written 
^(z)  - (W  + e F ^z)C(z)  = 0 


8ji?m*  dz^ 


(3.24) 


Fig.  3.5  Effect  of  electric  field  on  Cl-HHl  energy 
transition. 
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vrtiere  h is  the  Boltzmann  constant,  W is  the  energy,  and  is 

the  perpendicularly  applied  electric  field. 


The  zero-field  solutions  are  known,  i.e.. 


W = 


n 


nz 


(3.25) 


cos 


(3.26) 


It  is  therefore  natural  to  use  the  energy  of  the  ground  state 
as  the  unit  of  energy 


w. 


(3.27) 


and  as  unit  of  field  Fi  such  that  eFi  = Wi. 

In  order  to  transform  Equation  3.24  into  the  Airy 
differential  equation  it  is  convenient  to  use  a new  variable 

Z. 


Z = - 


8 31^ 


m 


(ehF,)' 


(W  + e F . z) 


(3.28) 


Note  that  W + e Fj^  z represents  the  kinetic 
8 71  m*/(ehFj^)  J 3 has  the  dimensions  of  energy. 


energy  and  that 
Z is  therefore 
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a normalized  kinetic  energy.  The  new  form  of  Equation  3.24 


can  be  written 

— - C(Z)  - Z C(Z)  = 0 
dZ 

Whose  solutions  are  of  the  form 

(3.29) 

^(Z)  = a Ai(Z)  + b Bi(Z) 

(3.30) 

Where  a and  b are  two  constants  and  Ai(Z)  and  Bi(Z)  are  the 


Airy  functions  [46].  The  boundary  conditions  are 
C(ZJ  = ^(Z_)  = 0 

(3.31) 

vdiere  Z^  = Z(±  ^ L^) . If  we  measure  the  energies  and  fields  in 


units  Wi  and  Fi,  i.e.. 

II 

(3.32a) 

Then 

(3.32b) 

Z^  = [-^]3(w±if) 

(3.33) 

and 


Ai(Z^)Bi  (Z  )-Ai(Z  )Bi(Z^)  = 0 
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(3.34) 


con5)letely  determines  the  eigen  energy.  Note  that  in 
the  normalized  form  the  effective  mass  and  the  well  thickness 
have  disappeared.  Therefore  Equations  3.33  and  3.34  are 
universal.  Figure  3.6  shows  the  universal  curve  for  the 
energy  for  the  first  state  in  an  infinite  quantum  well  in  the 
presence  of  a uniform  electric  field  perpendicular  to  the 
walls. 

The  InGaAs  wells  suffer  a conpressive  strain  vdien  grown 
on  GaAs  substrates  (AlGaAs  has  practically  the  same  lattice 
constant  as  GaAs)  due  to  the  fact  that  InGaAs  has  a higher 
value  of  the  lattice  constant.  The  presence  of  conpressive 
strain  and  the  fact  that  the  layers  are  grown  on  (lll)B  GaAs 
substrates  (as  was  previously  mentioned),  would  direct  the 
electric  field  vector  towards  the  substrate  (see  Fig.  3.7). 

Figure  3.8  shows  the  energy  band  structure  for  the  (100) 
and  the  (lll)B  grown  layers.  For  the  structures  that  we  are 
suggesting,  the  value  of  this  field  in  sample  SI  is 
calculated  to  be  1.1  10^  v/cm,  while  for  saitple  S2  it  is  5.5 
10^  V/cm.  Using  the  effective  infinite  well  method  described 
previously,  the  decrease  of  the  Cl-HHl  transition  energy  to  a 
lower  value  is  found  to  be  24  meV  and  11  meV  for  SI  and  S2, 
respectively. 


Normalized  Ground  State  Energu(NewvalueA)1(K'a1ue) 
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Fig.  3.6  Universal  curve  for  the  energy  of  the  first  state  in 
an  infinite  quantum  well  in  the  presence  of  a 
uniform  electric  field  perpendicular  to  the  walls  of 
the  well. 
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I Electric  Field 
f Inside  The 
Wells 


Fig.  3.7  Electric  field  direction  inside  the  wells. 


3.3.5  Conbined  Effects 


43 


Finally,  after  presenting  the  effects  of  the  effective 
mass,  the  strain,  and  the  presence  of  strain-induced  electric 
field  on  the  Cl-HHl  transition  energy,  we  add  these  effects 
together  for  the  (100)  and  the  (lll)B  grown  structures. 
Figure  3.9  shows  clearly  that  for  the  same  structures  grown 
on  (100)  and  (lll)B  GaAs  the  Cl-HHl  is  much  higher  in  the 
latter  due  to  higher  strain  induced  energy  shift. 

The  conputer  programs  written  to  calculate  the  Cl-HHl 
transition  energy  for  (100)  and  (111)  grown  MQWs  are 
presented  in  ^pendices  A and  B. 

3.3.6  Advantages  of  the  Proposed  Structures 

Most  of  the  reported  research  to  this  date  used 
InGaAs/GaAs  SQWs  [9,14,47]  with  thick  barriers,  eliminating 
strain  in  the  GaAs  barrier  layers.  We  propose  having  an 
InGaAs/AlGaAs  structure  instead  of  the  InGaAs/GaAs  structure 
because  in  the  latter  the  barriers  are  not  high  enough  to 
stop  the  tunneling  of  electrons  in  the  QWs  under  bias  and 
because  the  light  hole  energy  band  is  not  confined  within  the 
quantum  vrell  because  of  the  small  valence  band  offset  [48]. 

Furthermore,  the  use  of  AlGaAs  as  barriers  instead  of 
GaAs  assures  the  elimination  of  any  blurring  of  the  QW's 
absorption  by  the  absorption  of  the  barrier's  material  [49]. 

Under  bias,  the  Cl-HHl  transition  energy  moves  towards 


Energy 
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Growth 

Direction 


Fig.  3.8  (a)  Energy  band  structure  of  (100)  grown  samples. 

(b)  Energy  band  structure  of  (111)  grown  sanples. 
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Fig.  3.9  Combined  effects  determining  Cl-HHl  transition 
energy. 

(a)  For  structure  SI  (In  10%). 

(b)  For  Structure  S2  (In  5%). 
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higher  energy,  \^iile  the  barriers'  bandgap  transition  goes  to 
a lower  energy.  As  shown  in  Figure  3.10  the  GaAs  transition 
energy  can  cross  that  of  the  QWs  under  appropriate  bias, 
vdiile  for  AlGaAs  this  situation  is  not  likely  to  occur  until 
an  extremely  high  bias. 

In  the  next  chapters  we  will  investigate  these 
properties,  especially  the  presence  of  the  strain-induced 
electric  field. 
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Under  Bias  GaAs  and  AIGaAs  Trarsition 
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Fig  3.10  The  effect  of  electric  field  on  the  transition 

energy  in  the  viells'  in  S1-S2,  and  in  the  AlGaAs  and 
GaAs  layers. 


CHAPIER  4 

EFFECT  OF  OPTICAL  EXCITATION 
AND  MECHANICAL  STRESS 
ON  THE  CHARACTERISTICS  OF  THE  STRUCTURES 

4.1  Introduction 

One  of  the  methods  suggested  to  prove  the  existence  of  a 
strain-induced  electric  field  in  the  wells  is  to  use  PL 
measurement  of  the  QW  peak  position  and  its  shift  under 
increase  of  laser  intensity.  It  was  reported  [49]  that  the 
photocreation  of  electron-hole  pairs  would  cause  the 
screening  of  such  electric  field.  Moise  et  al.[49]  showed  a 
substantial  shift  of  the  strained  QW  peak,  grown  on  (lll)B 
GaAs,  to  a higher  energy  under  increased  optical  excitation. 

Tailoring  the  band  gap  transition  energies  can  be 
achieved  by  subjecting  a sample  to  mechanical  strain  vdiich 
would  play  the  same  role  as  the  one  played  by  the  growing 
mismatched  layers.  But  the  advantage  here  is  that  an  energy 
shift  can  be  controlled  by  placing  the  sanple  on  a transducer 
and  adjusting  the  applied  voltage  to  the  desired  value.  This 
controlled  transition  energy  shift  would  selectively  change 
the  band  gap  of  the  sample  thus  changing  its  optical 
characteristics.  The  increase  in  band  gap  due  to  ccxipressive 
strain  has  a reverse  effect  on  the  refractive  index.  While 
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under  tensile  strain  the  band  gap  decreases  in  um  increasing 
the  refractive  index.  We  will  investigate  the  effect  that 
applied  mechanical  tension  on  (111)  grown  sanples  has  on  the 
absorption  peak  position  and  on  its  intensity. 

4.2  Effect  of  Optical  Excitation 

Shining  a laser  beam  with  a photon  energy  higher  than 
the  expected  transition  energy  for  the  QWs,  will  cause  the 
creation  of  an  electron-hole  pair,  with  the  electron 
acquiring  enough  energy  to  be  transferred  to  the  conduction 
bcuid.  Once  there,  the  electron  is  bound  to  come  back  to  a 
lower  energy  state  passing  by  the  quantized  energy  states 
(back  to  the  valence  band)  releasing  in  the  process  a photon 
vdiose  energy  corresponds  to  the  energies  of  these  states  (see 
Fig.  4.1). 

The  photons  created  in  the  process  are  passed  through  a 
monochrometer  that  selects  the  appropriate  light  wavelength, 
via  a grating.  These  photons  having  a certain  wavelength  are 
then  detected  by  a thermoelectrically  cooled  photomultiplier 
tube.  The  latter  is  connected  to  a computer  where  the 
information  about  the  intensity  and  the  wavelength  of  the 
radiation  is  stored. 

As  the  intensity  of  the  applied  laser  light  is 
increased,  the  number  of  created  electron-hole  pairs 
increases  given  rise  to  an  electric  field  that  would 
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partially  screen  the  strain- induced  electric  field  present  in 
the  QWs.  As  was  mentioned  in  the  previous  chapter,  the 
presence  of  an  electric  field  inside  a QW  will  cause  the  Cl- 
HHl  transition  energy  to  be  reduced.  So  an  increase  in  the 
intensity  of  the  optical  excitation  would  yield  a shift  of 
the  PL  peak  position  of  the  QW  to  a higher  energy  only  in 
sanples  vdiere  a built-in  electric  field  (111  grown  samples) 
already  exists. 

Moise  et  al.  [49]  reported  that  a strain  induced 
electric  field  in  the  QWs  of  9x10^  V/cm  corresponds  to  an 
interface  charge  density  of  6.6xl0ii  cnr^.  Therefore,  for  the 
structures  at  hand,  the  density  of  carriers  needed  for 
complete  screening  of  the  internal  fields  comes  out  to  be 
8.067x1011  cin-2  for  the  sample  SI  (10%  In),  and  4.034xl0H  cm-2 
for  sairple  S2  (5%  In). 


4.3  PL  Measurements 


The  Photoluminescence  setup  used  in  the  measurements  is 
depicted  in  Fig.  4.2.  The  PL  was  measured  at  13  K using  a 
0.64  m spectrometer.  A 40  mw  Ar‘*'-ion  laser  operating  at 
514.5  nm  was  used  as  the  excitation  source  and  the  signal  was 
detected  with  a thermoelectrically  cooled  photomultiplier. 
All  the  measurements  were  done  under  identical  conditions  for 
all  sanples  under  study  in  order  to  conpare  the  results. 

From  the  results  presentedin  Figs. 4. 3-6,  it  can  be  seen 
that  in  the  case  of  the  (111)  grown  layers  the  intensity  of 
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Fig.  4.1  The  e~  transition  frcxn  Cl  to  HHl. 
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the  MOW  luminescenoe  is  1.2  to  2 times  stronger  than  that  of 
the  (100)  grown  layers,  and  the  FWHM  for  (111)  is  smaller 
than  that  of  the  (100)  grown  MQWs,  implying  an  increase  in 
the  optical  matrix  element  [7,14].  This  shows  that  the 
optical  quality  of  (111)  grown  material  is  excellent  and  is 
better  than  the  layers  grown  on  (100)  [18]. 

when  increasing  the  excitation  power  the  number  of 
carriers  inside  the  wells  increases  causing  the  screening  of 
the  strain-induced  electric  field.  This  will  cause  the 
energy  bands  to  move  closer  to  a flat  band  case  producing 
higher  luminescence  intensity  and  lower  FWHM.  The  PL 
measurements  depicted  in  Figs  4.3,5  show  a shift  of  the  MQW 
peak  for  the  (111)  samples  toward  higher  energies  and  a 
decrease  in  the  FWHM  with  increasing  laser  intensity  while 
the  same  trend  is  absent  for  (100)  (Figs.  4.4,6).  This 
confirms  the  presence  of  a built-in  electric  field  inside  the 
(111)  grown  layers. 

Using  the  Kronig-Penney  model  presented  in  Chapter  3,  we 
generated  Figs.  4.7  and  4.8  that  show  the  change  in  the  MQW 
peak  energy  when  the  strain-induced  electric  field  is 
screened  for  the  structures  SI  and  S2,  respectively.  For  the 
levels  of  power  used  in  the  measurements,  the  shifts  of  MQW 
are  measured  to  be  2.44  meV  and  0.5  meV  for  SI  and  S2, 
respectively. 


4.4  Amount  of  E-Field  Screening 
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In  order  to  determine  the  intensity  of  the  light 
reaching  the  wells,  we  need  to  know  the  laser  intensity  at 
the  surface  of  the  sample,  the  power  reflection  from  both 
sides  of  the  glass  window  of  the  cryostat  and  from  the  top 
of  the  san^le  itself,  and  the  absorption  of  the  layers  atop 
of  the  wells'  region. 

The  equation  for  the  reflectivity  at  normal  incidence  is 
given  by  [50], 


where  ni  and  n2  are  the  refractive  indices  of  the  adjacent 
regions  (air-glass,  glass-air,  and  air-sanple) . 

For  the  cryostat  glass  window,  the  refractive  index  is 
1.5  yielding  a reflectivity  of  0.04  for  each  of  the  sides  of 
the  glass.  The  refractive  index  of  GaAs  is  4.14  (see  Chapter 
5 for  refractive  index  calculations)  at  514.5  nm  producing  a 
reflectivity  of  0.373.  The  ccxnbination  of  all  these  factors 
will  reduce  the  intensity  by  a factor  of  0.5778  upon  reaching 
the  sanple's  surface. 

To  be  able  to  change  the  intensity  of  the  laser 
excitation  neutral  density  filters  were  placed  in  the  path  of 
the  laser  light.  Full  power  in  the  case  of  SI  is  24.58  mw, 
with  the  neutral  density  filter  NDl  it  is  2.458  nW  and  with 
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Fig  4.2  PL  setup 
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ND2  it  is  0.5  mW.  For  S2,  Full  povrer  is  13.76  mW,  NDl  is 
1.87  and  ND2  is  0.5  ntN. 

The  beam  coming  out  of  the  laser  is  considered  to  be  a 
Gaussian  beam  with  a TEMq,o  mode  [51].  The  equations 
governing  the  spot  size  for  such  a beam  are  : 


w2(z)  = 


(4.2) 


z 


0 


njtw^ 

O 

X 


(4.3) 


vdiere  Wq  is  the  spot  radius  at  a position  z=0,  w(z)  the  spot 
size  at  a position  z,  and  n is  the  medium's  refractive  index, 
and  Wg  at  the  sanple's  surface  was  determined  to  be  36  jm. 
Using  the  above  information,  the  laser's  intensity  can  be 
derived  from  : 

P 

T = 8_  , . 


where  Ig,  Pg  and  Wg  are  the  laser  intensity,  power  and  spot 
size  at  the  surface  of  the  sanple,  respectively. 

After  combining  these  different  factors,  we  can 
calculate  the  intensity  of  the  light  hitting  the  sample's 
surface.  For  full  laser  power  in  the  case  of  SI  the 
intensity  is  571.51  W cm-2,  with  NDl  in  place  it  is  57.15  w 
cm-2,  and  with  ND2  in  place  it  is  11.62  w cm-2.  m the  case 
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of  S2  Full  corresponds  to  an  intensity  of  320  W cm-2,  with  NDl 
43.48  W cnr2,  and  with  ND2  11.62  W cnr^. 

The  absorption  loss  caused  by  the  presence  of  layers 
atop  of  the  C2W  region  is  given  by  [52] 

I = lo  e-oL  (4.5) 

1 

a = A (hv  - Eg)  2 (4.6) 

vdiere  I is  the  final  intensity,  Iq  is  the  initial  intensity, 
a is  the  absorption  coefficient,  A is  a constant,  hv  is  the 
energy  of  the  laser  light.  Eg  is  the  bandgap  energy,  and  L is 
the  width  of  the  layer. 

We  assume  that  all  the  e-h  pairs  created  in  the  p+ 
layers  will  be  collected  in  the  i region.  Considering  that 
the  spot  size  does  not  change  significantly  vdien  traveling 
through  these  layers,  the  intensity  reaching  the  wells 
suffers  some  absorption  in  the  AlGaAs  layer  atop  the  QW 
region  (reduction  of  40%).  For  sample  SI  the  intensity 
reaching  the  wells  is  about  157  w citr^  and  for  sanple  S2  it  is 
about  163  W cm-2  for  Full  power. 

Moise  et  al.  [49]  showed  that  at  a wavelength  of  514.5 
nm  the  intensity  needed  to  create  2x1 0^2  cm-2  e-h  pairs  is 
413.6  W cm-2.  In  our  case,  200  (571.51  x 0.5778  x 0.6)  W cm-2 
would  yield  9.65xl0H  cm-2  e-h  pairs,  and  112  (320  x 0.5778  x 
0.6)  W cm  2 would  yield  5.41  lOHxcm-2  e-h  pairs,  causing  a 
total  screening  of  the  electric  field. 


Photon  Countt/s 


Fig  4.3  C3W  peak  shift  with  increased  laser  power. 


Fig  4.4  QW  peak  shift  with  increased  laser  power 
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Fig  4.5  QW  peak  shift  with  increased  laser  power. 
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Fig  4.6  QW  peak  shift  with  increased  laser  power 
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Fig.  4.7  Effect  of  screening  the  E field  in  the  C2WS. 


Electric  Field  Present  in  MQWs  ( 1 oV/cm) 


Fig.  4.8  Effect  of  screening  the  E field  in  the  QSfe 


4.5  Mechanical  s-hrain  Effect 
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In  order  to  investigate  the  effect  that  applied 
mechanical  strain  has  on  the  optical  characteristics  of  (111) 
grown  samples,  a system  of  measurement  was  set  up  (Fig.  4.9). 
The  set  up  was  placed  inside  the  cryostat.  The  hole  going 
through  the  setup  is  aligned  with  a hole  on  the  tip  \diere  it 
is  placed.  A vdiite  light  source  is  placed  in  the  back  of  the 
cryostat  with  a lens  to  focus  this  light  into  the  desired 
hole.  The  light  transmitted  through  the  sample  is  then 
collected  by  the  monochrcxnBter-PMT  system,  in  order  to  obtain 
the  position  of  the  Cl-HHl  absorption  peak. 

The  samples  were  placed  on  a diaphragm  which  can  be 
moved  forward  by  the  force  of  an  electromagnet  applying  a 
tensile  strain  on  the  sanples  (Fig.  4.10).  Increasing  the 
applied  voltage  in  turn  causes  an  increase  in  the  applied 
tensile  strain. 

The  structure  under  investigation  was  grown  on  (lll)B 
GaAs  and  is  shown  in  Fig.  4.11.  it  has  20  quantum  wells  made 
of  lno.12Gao.88As/  Alo.15Gao.85As  with  a well  width  of  110  A, 
and  a barrier  width  is  200  A. 

The  effect  of  increasing  the  coil  voltage  on  the  Cl-HHl 
transition  for  the  sanple  at  room  teirperature  and  at  13  K is 
shown  in  Figs.  4.12-13.  It  can  be  seen  that  this  transition 
moves  towards  a lower  energy  suggesting  the  presence  of  a 
tensile  strain  affecting  the  sanple.  This  decrease  in  energy 


is  accompanied  by  a reduction  in  the  absorption  peak 
intensity. 

The  trend  is  the  same  for  both  ten^jeratures  but  for  the 
measurement  done  at  13  K the  energy  shift  is  larger  than  that 
at  room  ten^rature  (16.4  meV  at  300  K,  47.5  meV  at  13  K). 
This  Ccin  be  explained  by  the  fact  that  the  adhesive  used  to 
place  the  sample  on  the  diaphragm  (Crazy  Glue),  shrinks  at 
low  temperatures,  adding  more  tension  to  the  sanple,  moving 
the  Cl-HHl  transition  higher  in  energy. 

Care  was  giving  to  the  reduction  of  the  area  vrtiere  the 
adhesive  material  makes  contact  with  the  san^le.  It  was 
found  that,  for  13  K measurement,  having  a large  contact  area 
will  ccxipletely  bleach  the  Cl-HHl  absorption  peak. 

To  learn  more  about  this  applied  tensile  stress,  a 
simulation  was  run  using  the  modified  Kronig-Penney  model 
(including  electric  field)  introduced  in  Chap  3,  to  evaluate 
the  amount  of  conpressive  strain  and  electric  field  remaining 
in  the  layers.  Figs.  4.14  and  4.15  show  the  amount  of  the 
expected  change  in  the  lattice  constant  and  the  strain- 
induced  electric  field.  Increasing  the  coil  voltage  will 
reduce  the  percentage  lattice  mismatch  and  the  electric  field 
These  values  are  different  for  the  room  tenperature  and  13  K 
cases  because  of  the  reasons  mentioned  previously. 

The  application  of  a coil  voltage  of  2 volts  will  reduce 
the  strain  by  0.40%  and  0.51%,  and  the  electric  field  by 
6.195x104  V cm-l  and  8.08x104  v cm-l,  in  the  cases  of  room 
tenperature  and  13  K.  That  is  a 47.66%  and  46.3%  reduction 


62 


Strained 


MQW 

(111) 


Diaphragm 


Supporting 

Post 


Variable 
Power 
Supply 
for  Coil 


Yoke 


Coil 


Fig.  4.9  The  set  up  used  to  apply  machanical  strain. 
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in  the  strain  and  the  electric  field,  respectively,  for  the 
room  ten^rature,  \diile  the  reduction  is  73.38%  and  70.87%, 
respectively,  for  the  13  K case. 

In  conclusion,  vre  have  proved  the  presence  of  a strain 
induced  electric  field  in  the  (111)  grown  sanples  by  the  use 
of  intensity  controlled  photoluminescence.  These  same 
measurements  showed  no  such  electric  field  in  (100)  grown 
samples.  Then  the  screening  of  this  electric  field  was 
determined  from  these  previous  measurements. 

Finally,  the  application  of  mechanical  strain  on  a (111) 
grown  saiq)le  showed  a substantial  change  in  the  position  of 
Cl-HHl  absorption  peak  under  low  bias  values.  This  feature 
can  be  used  to  design  and  operate  (111)  based  switches. 
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Fig  4.10  Sanple  under  strain. 
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Fig  4.11  Structure  of  sample  for  strain  studies. 
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Fig  4.12  Room  tenperature  Cl-HHl  absorption  peak  shift  under 
applied  strain. 
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Fig  4.13  13  K Cl-HHl  absorption  peak  shift  under 
applied  strain. 


1.6 


2 


Luminescence  Peak  (eV) 


68 


03 

03 

0.7 

0.6 

03 

0.4 

03 


Fig.  4.14  Simulation  of  effect  of  applied  strain  on  the 
built-in  electric  field  and  on  the  strain 
percentage  in  the  layers  at  room  tenperature. 
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Fig.  4.15  Simulation  of  effecft  of  applied  strain  on  the 
built-in  electric  field  and  on  the  strain 
percentage  in  the  layers  at  13  K. 
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CHAPTER  5 

CHARACTERIZATION  OF  IHE  THIRD 
ORDER  NONLINEARITY 


5.1  Intrcxiuction 


Intensity  dependent  nonlinearities  have  been 
investigated  for  AlGaAs/GaAs  and  strained  InGaAs/AlGaAs  MQW 
structures  grown  on  (100)  GaAs  substrates.  Different 
experiments  can  be  performed  to  determine  the  nonlinear 
refractive  index  ti2,  using  the  pun^)  and  probe  method  or  the 
hjiarid  Mach-Zender  interferometer  method  [53-55]. 

In  this  chapter  we  present  the  values  obtained  for  the 
nonlinear  refractive  index  n2  in  the  structures  grown  on 
(100)  and  (lll)B  GaAs  substrates.  A hybrid  Mach-Zender 
interferometer  was  used  to  experimentally  determine  the  value 
of  n2 . This  work  was  limited  to  TE  CW  measurements  at 
wavelengths  of  980  nm  and  992.2  nm. 

The  structures  (S1,S2)  were  designed  to  insure  a single 
mode  slab  waveguide  operation  at  the  above  mentioned 
wavelengths,  where  the  guided  mode  peak  has  a large  overlap 
with  the  10  period  MQW  region. 
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5.2  Wavecfiiidina  Properties  of  the  Structures 
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5.2.1  Multilayer  Matrix  Theory 

In  order  to  determine  the  appropriate  conpositions  of 
the  different  layers  to  insure  a single  mode  operation  the 
multilayer  matrix  method  was  used  [56,57].  This  method 
depends  on  equating  the  fields  at  the  boundaries  of  each 
layer  and  deriving  the  dispersion  relation  for  the  total 
structure. 

The  theory  of  multilayer  stacks  [56,57]  defines  two 
field  variables  U and  V given  by 

U = Ey  (5.1a) 

V = (5.1b) 

which  describe  the  transerve  variation  of  the  optical  field. 

U and  V are  the  transverse  field  distribution  in  a particular 
layer  of  constant  refractive  index  n.  These  definitions  are 
chosen  because  U(x)  and  V(x)  are  quantities  that  are 
continuous  at  the  layer  boundaries.  From  the  following 
equations,  derived  from  Maxwell's  equations 

pEy  = -<omH^ 


(5.2) 
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dE 


dx 


- = -j(«nH 


(5.3) 


dx 


+ jPH*=  -jw£E, 


(5.4) 


we  certain  the  relations 


U = -j  V 


v‘  = j(p^-  k^)U 


(5.5a) 

(5.5b) 


where  the  prime  indicates  differentiation  with  respect  to  x. 
Both  u and  V c*)ey  the  transverse  wave  equation 


ti  = (p^-  n^  k^)u 


(5.6) 


The  genercil  solution  of  the  wave  equation  in  this  layer  is  : 


U = A exp(-jKX)  + B exp(jKX) 


V = K [a  exp(-jKx)  - B exp(jKX)  ] 


(5.7a) 

(5.7b) 


where 


= n^k^  - p^ 


(5.8) 
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The  A eind  B constants  can  be  determined  from  the  values  Uq  = 
U(0)  and  Vo  = V(0)  at  the  input  plane  x = 0 of  the  layer,  so 
vie  have 


A = L (u  + 

(5.9) 

B = 1 (u,  - 

2 ' ° K / 

(5.10) 

A rearrangement  of  equations  (5.7-10)  leads  to  a sinple  matrix 
relation  between  the  output  quantities  U,  V and  the  input 
quantities  Uq  and  Vq 


where  the  pairs  (Uo,Vq)  and  (U,V)  have  been  written  as 
vectors,  and  M is  the  characteristic  matrix  of  the  layer. 
This  matrix  has  the  form 


cos(kx) 

(i)sin(KX) 

jKsin(KX) 

COS(KX) 

(5.12) 


Next,  consider  a stack  of  n layers  sandwiched  between 
substrate  and  cover  as  seen  in  Fig.  5.1.  The  layers  are 
labeled  starting  from  the  substrate,  the  layer  thicknesses  are 
hi  and  the  layer  indices  are  ni,  vdiere  i = 1 to  n.  The  output 
field  variables  for  each  layer  are  Ui  and  Vi. 
characteristic  matrices  at  the  layers  are 


The 
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cos(K .h^) 

(i)sin(Kihi) 

M.  = 

1 

jK^sin(K^h^) 

cos(K  ^h^) 

vdiere 


(5.13) 


(5.14) 


The  corresponding  field  variables  are  related  by 


(Ui-i] 

Wi-i) 


(5.15) 


Using  matrix  multiplication,  a simple  relation  between  the 
input  variables  Uq,  Vq  at  the  substrate  and  the  output 
variables  Un,  Vn  at  the  cover  can  be  obtained 


where  M is  the  characteristic  matrix  of  the  stack.  It  is 
given  by  the  product  of  the  individual  layer  matrices 


M 


®21 


® 12 
®22 


= M^.M2. 


(5.17) 


vdiere  mu,  mi2,  m2i  and  11^2  are  the  matrix  elements. 
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The  reflection  and  transmission  coefficients  for  the 


light  incident  on  the  stack  can  be  determined  as  follows, 

Assuming  the  light  is  incident  fran  the  substrate  side,  it  can 
described  by  Eqs.  (5.7a)  and  (5.7b)  with 

P'  = -Ys  (5.18) 

along  with  the  incident  amplitude  Ag  and  the  reflected 
anplitude  Bg.  Therefore,  the  ir^ut  variables  are 

U„  = + B,  (5.19a) 

''o  = ■'J*.  - (5-19b) 


On  the  cover  side,  the  transmitted  light  has  an  anplitude  Ac, 
and  the  corresponding  output  variables 

Un  = K (5.20a) 

Vn  = (5.20b) 


vdiere 


(5.21) 


The  optical  field  in  the  substrate  and  the  cover  layers 
consist  of  evanescent  waves.  The  corresponding  fields  are 
written  in  the  form 
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z 


Fig.  5.1  Schematic  of  n layer  stack. 


U = A exp(YX)  + B exp(-YX) 
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(5.22) 

V = jY(A  exp(YX)  - B exp(^x)  ) 

(5.23) 

which  obey  the  transverse  wave  Eqs.  (5.6).  Instead  of  the 
transverse  propagation  constants  and  the  transverse 
decay  constants  Vg  Yc  should  be  used.  Considering  that 
the  fields  should  decay  away  from  the  multilayer,  the  input 


and  output  field  variables  are 

(5.24a) 

'^0  = iVs*. 

(5.24b) 

(5.24c) 

'^n  = -3^0*0 

Inserting  Eqs.  (5.24a-d)  into  equation  5.16  yields 

(5.24d) 

As  = (“^ii  - :)Yc"‘i2)Bc 

(5.25) 

jYgAg  = (m2i  - jYcm22)Bc 

(5.26) 

Dividing  these  two  relations  leads  to  the  dispersion  relation 
for  a multilayer  slab  guide 


3(Ygmii  + Yc"^22)  = "^21  - YgYcroi2 


(5.27) 
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This  relation  establishes  the  connection  between  the  frequency 
(0  of  the  light  and  the  propagation  constant  p of  the  mode 

guided  by  the  multilayer. 

5.2.2  Refractive  Index  of  Different  Layers 

The  refractive  index  value  for  GaAs,  AlGciAs  and  InGaAs 
have  been  reported  earlier  [58-60].  These  values,  at  a given 
wavelength,  can  be  derived  by  the  use  of  the  single-effective- 
oscillator  (SEO)  method  given  by  Nahory  et  al.  [60],  Utaka  et 
al.  [61],  and  Brcberger  and  Lingren  [62] 


The  parameters  Eq,  Ed  and  Eg  are  expressed  in  ev  by  the 
following  enpirical  relations  : 


n = 


V 


(5.28) 


Where 


(5.29) 


i)  For  AlxGai_xAs 
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Eo  = 3.65  + 0.871  X + 0.179  x2 

(5.30) 

Ed  — 36.1  - 2.45  X 

(5.31) 

Eg  = 1.424  + 1.266  x + 0.26  x2 

(5.32) 

Eo  = 3.65  - 2.15  y 

(5.33) 

Ed  = 36.1  - 19.9  y 

(5.34) 

Eg  = 0.36  + 0.505  (1-y)  + 0.555  (l-y)2 

(5.35) 

Replacing  Eqs.  (5.30-32)  in  equation  will  yield  a value 
of  the  refractive  index  with  respect  to  energy  for  different 
values  of  the  A1  conposition  x. 

For  InGaAs  the  shift  in  energy  due  to  lattice  mismatch 
growth  would  change  the  refractive  index  value  considerably 
for  high  Indium  percentage.  Das  and  Bhattacharya  [63]  has 
showed  a decrease  in  the  refractive  index  of  InGaAs /GaAs 
heterostructures  due  to  the  increase  in  the  bandgap  of  InGaAs 
material  under  conpressive  strain.  For  our  sanples  the  value 
of  the  Indium  percentage  (10%  in  SI,  5%  in  S2)  is  too  small  so 
this  effect  is  not  present.  The  refractive  index  of  GaAs, 
Alo.2Gao.8As,  and  lno.05Gao.95As  are  presented  in  Fig.  5.2. 
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The  refractive  index  of  the  MQW,  given  by  Glick  and 
Reinhart  [64],  is  derived  from 


Using  the  determined  values  of  the  refractive  indices  in 
the  equations  derived  from  the  multilayer  matrix  theory 
described  above,  the  modefield  can  be  obtained.  These 
modefields  are  shown  in  Figs.  5.3  to  5.6  for  the  TEq  mode.  In 
these  figures  the  position  zero  corresponds  to  the  bottom  of 
the  QW  region.  It  can  be  seen  that  the  field's  peak  is  within 
the  MQW  region. 


5.3.1  Introduction 

A linear  dielectric  medium  is  characterized  by  a linear 

relation  between  the  polarization  density  and  the  electric 
field,  P = E , where  £q  is  the  permitivity  of  free  space 

and  X is  the  electric  susceptibility  of  the  medium.  A 
nonlinear  dielectric  medium,  on  the  other  hand,  is 
characterized  by  a nonlinear  relation  between  P and  E. 

Since  externally  applied  optical  electric  fields  are 
typically  small  in  conparison  with  characteristic  interatomic 
or  crystalline  fields,  the  nonlinear  term  is  usually  weak. 


(5.36) 


5.3  NOnlinearities  in  Proposed  Structures 
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The  relation  between  P and  E is  then  approximately  linecir  for 
small  E,  deviating  only  slightly  from  linearity  as  E 
increases.  Under  these  circnimstances , it  is  possible  to  expand 
the  function  that  relates  P to  E in  a Taylor's  series  about  E 
= 0, 

P = a,E  + 1-  a,E^  + ^ a,E^  + ...  (5.37) 

The  coefficients  al,  a2,  and  a3  are  the  first,  second,  and 
third  derivatives  of  P with  respect  to  E at  E = 0.  The  first 
term,  vdiich  is  linear,  dominates  at  small  E.  The  second  term 
represents  a quadratic  or  second-order  nonlinearity,  the  third 
term  represents  a third  order  nonlinearity,  and  so  on. 

The  above  equation  is  usually  written  as  [65] 


In  media  possessing  centrosymmetry,  the  second-order  term  is 
absent  since  the  polarization  must  reverse  exactly  when  the 
electric  field  is  reversed.  The  dcxninant  nonlinearity  is  then 
of  third  order. 


P = EpX  E + 2d  E^  + E^  + ... 


(5.38) 


(5.39) 


A material  with  large  third  order  nonlinearity  is  sometimes 
called  a Kerr  medium. 
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u 


Fig.  5.2.  The  refractive  index  vs.  Energy  for  GaAs,  Alo.2Gao.9As,  and  Ino.05Gao.95As 
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Fig.  5.3.  The  modefield  for  the  structure  SI  at  = 980 


Mode  Field  in  SI  at  0.9922 
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Fig.  5.4.  The  modefield  for  the  structure  SI  at  Xo  = 992.2 


Mode  Field  in  S2  at  0.98  /<m 
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Fig.  5.5.  The  modefield  for  the  structure  S2  at  Xo  = 980 
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Fig.  5.6.  The  itodefield  for  the  structure  S2  at  Xo  = 992.2 
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In  accordance  with  Eq.  5.39,  the  response  of  a third-order 
nonlinear  medium  to  a monochromatic  optical  field  E(t)  = Re 
{E(  (u)exp( j cot) } is  nonlinear  polarization  PuL(t)  containing  a 
cxanponent  at  frecjuency  co  and  another  at  frecjuency  3 co, 

= 3 I E(co)  I ^ E(co)  (5.40) 


Pnl(3«)  = x‘^^E^(o) 


(5.41) 


The  polarization  component  at  frequency  co  in  Eq.  5.40 
corresponds  to  an  incremental  change  of  the  susceptibility  Ax 

at  frecjuency  co  given  by 


= 

where 


E(co) 


j ^ I E(co) 


= 3 X 
2 


2t) 

wave.  Since  n^  = 1 
refractive  index  An 


(3) 

is 

+ 


I E(co)  1^  = 6 x^^’tiI  (5.42) 

the  optical  intensity  of  the  initial 

X,  this  is  ecjuivalent  to  an  incremental 

, Ax 

— Ax  = — , so  that 

\^l  2n 


An=  ^ 


X^^^I  = n 


(5.43) 


Thus  the  third  order  nonlinearity  causes  a change  in  the 
refractive  index  proportional  to  the  optical  intensity.  The 
overall  refractive  index  is  therefore  a linear  function  of  the 
optical  intensity  I, 


n(I)  = n + n2l 
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(5.44) 


(5.45) 


is  called  the  nonlinear  refractive  index. 

This  effect  is  known  as  the  optical  Kerr  effect.  This  effect 
is  self-induced  in  that  the  phase  velocity  of  the  wave  depends 
on  tte  wave's  own  intensity. 

5.3.2  Measurements  of  the  Nonlinear  Refractive  Index 

Consider  a waveguide  vdiere  an  optical  field  is  launched. 
Assuming  a linear  dependence  of  the  phase  shift  on  the  input 
povrer  [55],  one  can  express  the  value  of  n2  as 


vdiere  A(j)  is  the  relative  phase  shift.  Pi  is  the  input  power, 
A is  the  guide  cross  section,  Iz  is  the  length  of  the 


light,  and  ct  is  the  loss  coefficient  of  the  waveguide. 
Thues,  by  measuring  the  phase  shift  A(|>  and  the  loss 
cosfficient  a,  the  value  of  U2  can  be  determined. 


(5.46) 


waveguide,  Xo  is  the  free-space  wavelength  of  the  incident 
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The  Mach-Zender  interferometer  is  shown  in  Fig.  5.7, 
with  the  2mm  long  waveguides  placed  in  one  of  the  cirms.  The 
interference  fringes  produced  due  to  the  difference  in  the 
optical  path  length  between  the  two  arms  of  the  Mach-Zender 
are  captured  on  a computer  using  a CCD  camera  and  a Frame 
Grabber  Card.  Increasing  the  input  power  causes  a shift  in 
the  interference  pattern  (Fig.  5.8).  The  peeUcs  in  Fig.  5.8 
correspond  to  the  bright  areas  in  the  interference  pattern, 
while  the  valleys  correspond  to  the  dark  areas.  The  phase 
difference  corresponding  to  two  consecutive  peaks  (valleys) 
is  180  degrees.  From  that  we  can  calculate  the  phase  shift 
A(j)  with  power  vdien  displaying  the  interference  pattern  for 
two  power  values.  The  phase  shift  with  the  increase  of  input 
power  is  presented  in  Figs. 5. 9 and  5.10  for  all  four  sanples 
at  Xo  = 980  and  992.2  nm.  It  can  be  seen  that  SI  (111) 
sanple  shows  a leveling  off  at  high  optical  powers  suggesting 
a high  optical  powers  suggesting  a saturation  of  the 
nonlinearity. 

The  value  of  a for  each  of  the  sanples  was  obtained  by 
the  cut  back  method.  The  estimated  values  of  TI2  for 
different  powers,  for  wavelengths  of  980  and  992.2  nm,  for 
all  the  san^jles  are  given  in  Table  5.1  and  5.2.  The  values 
of  n2  obtained  for  the  (100)  grown  sanples  are  ccwiparable  to 
those  reported  in  refs.  [9,66,67].  These  data  show  that  the 
(111)  grown  sanples  have  30-50%  higher  nonlinear  refractive 
index  conpared  to  their  (100)  grown  counterparts. 
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The  origin  of  112  can  be  thermal  or  electronic.  For  CW 
measurements  both  of  them  are  present,  while  in  pulsed 
measurements,  with  pulse  duration  less  than  10  ns,  the  origin 
of  n2  is  mainly  electronic.  The  thermal  effect  moves  the 
transition  energies  of  the  MQWs  to  lower  values  increasing 
the  sanples'  refractive  index  (n2  is  positive).  On  the  other 
hand,  the  electronic  effect  will  yield  an  increase  in  the 
value  of  the  transition  energies,  decreasing  the  refractive 
index  (n2  is  negative). 

The  transition  energy  for  SI  (100)  is  from  Fig.  (3.9)  is 
1.3716  eV,  while  for  Sl(lll)  it  is  1.3974  eV.  For  S2(100) 
the  value  is  1.4019  eV,  and  for  S2(lll)  it  is  1.4298  eV.  Bt»r 
Xo  = 980  nm  (1.2653  eV),  the  energy  difference  between  the 
laser  and  the  closest  Cl-HHl  transition  is  106.29  meV 
(farthest  is  164.5  meV),  vrtiile  for  Xo  = 992.2  nm  (1.2497  eV) 
the  closest  becomes  121.85  meV,  while  the  farthest  will  be 
180  meV.  For  sources  operating  so  far  from  the  transition 
energies,  the  thermal  effect  is  not  large,  owing  the  origin 
of  the  nonlinearity  to  the  electronic  effect (n2  is  negative). 

For  (100)  saitples  only  the  electronic  effect  is  present, 
while  for  the  (111)  grown  samples  there  is  an  additional 
factor  influencing  the  value  of  n2.  That  factor  is  the 

screening  of  the  strain-induced  electric  field  increasing  the 
value  of  n2. 

The  value  of  n2  is  closely  connected  to  that  of  a, 
increasing  a would  increase  n2.  Eventhough  the  (100)  sanples 
are  closer  in  energy  to  the  energy  of  the  lasers  used  in  this 
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experiment,  the  value  of  a is  higher  in  the  (111)  samples. 


This  is  probably  due  to  the  dispersion  along  the  layers 
interfaces  in  these  sanples. 

The  values  of  n2,  for  (100)  san^les,  reported  here  are 
conparable  to  values  given  by  others  [63].  The  decrease  in 
the  value  of  n2  between  Xo  = 980  and  992.2  nm  can  explained 
by  the  reduction  of  the  absorption  coefficient  vdien  operating 
farther  away  from  the  band  edge.  This  trend  was  observed  by 
other  researchers  doing  wavelength  dependent  measurements  of 
n2[68]. 

Furthermore,  Moise  et  al.  [67]  using  time-resolved 
photoluminescence  measurements  of  similar  p-i-n  structures 
grown  on  (100)  and  (111)  substrates  with  a single 
InGaAs/AlGaAs  quantum  well  in  the  i region,  showed  that  the 
(111)  grown  sample  exhibited  a smaller  luminescence  decay 
time  than  the  (100)  grown  sample  in  response  to  an  increase 
in  photogenerated  carrier  density  (1.625  times  shorter). 
This  shows  that  future  pulsed  measurements  may  demonstrate 
much  better  performance  of  (111)  strained  MQWs  nonlinearities 
than  that  of  the  (100)  grown  MQWs. 

The  figure  of  merit  for  the  (100)  and  (111)  sanples  is 
given  by 


TO 


vAiere  x is  the  decay  time, 
so. 


(5.47) 
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F 


(111)  _ ^2(111)'*^  (100)“  (100) 
(100)  "2(100)''(111)“(111) 


F 


(5.48) 


for  SI 


= 1.69 


(5.49) 


and  for  S2 


= 0.926 


(5.50) 


The  70%  inprovement  in  the  figure  of  merit  value  in  the 
case  of  SI  (10%)  makes  the  (lll)B  grown  strciined  MQWs  with  a 
high  In  percentage  a better  choice  for  the  design  and 
fabrication  of  nonlinear  switches.  While  for  S2(5%  In)  the 
figure  seems  to  be  the  same,  showing  that  the  screening 
effect  here  does  not  play  as  a bigger  role. 
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Fig  5.7.  The  Mach-Zender  interferaneter  setup. 
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U) 


Fig.  5.8.  The  shift  in  interference  pattern  with  increase  of  input  laser 
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Fig.  5.9.  Relative  phase  shift  vd.th  increase  in  input  power 
for  Xo  = 980  nm. 
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Fig.  5.10.  Relative  phase  shift  with  increase  in  input  power 
for  Xo  = 992.2  ran. 


Table  5.1  n2  for  X©  = 980  nm 
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TE 

X = 0.980 

Length  (cm) 

Loss 

Coefficient  a 
(cirri) 

n2  (cm^A?) 

Sl(lll) 

0.200 

7.3 

8.69  10-8 

Sl(lOO) 

0.225 

5.7 

6.55  10-8 

82(111) 

0.200 

4.3 

4.83  10-8 

82(100) 

0.250 

1.5 

3.05  10-8 

Table  5.2  xi2  for  Xo  = 992.2  nm 


TE 

X = 0.9922 

Length  (cm) 

Loss 

Coefficient  a 
(cm-i) 

n2  (cm2A^) 

81(111) 

0.200 

1.5 

3.52  10-8 

81(100) 

0.225 

1.2 

2.65  10-8 

82(111) 

0.200 

0.9 

2.64  10-8 

82(100) 

0.250 

0.3 

1.70  0-8 

CHAPTER  6 

COJCLUSIONS  AND  RECOMMENDATIONS 

The  Study  of  the  proposed  structures  grown  on  (lll)B 
GaAs  yielded  some  inportant  features.  The  first  of  vdiich  is 
the  proof  of  the  existence  of  a strain-induced  electric  field 
in  the  strained  wells  using  photoluminescence  measurement. 
This  comes  from  the  fact  that  increasing  the  optical 
excitation  intensity  used  to  test  the  sanples  would  increase 
the  e-h  pairs  created  in  the  wells.  Since  the  presence  of  an 
electric  field  in  the  wells  would  cause  a red  shift  (QCSE)  in 
the  transition  energies,  the  e-h  pairs  would  act  to  screen 
any  existing  built-in  electric  field  reducing  the  effect  that 
such  field  has  on  the  transitions  energies.  That  means, 
under  increased  excitation  the  transitions  energies  should 
move  to  higher  values.  That  is  what  was  observed  in  the  PL 
spectra  for  (lll)B  grown  samples,  vdiile  this  trend  was  absent 
in  the  (100)  grown  samples. 

The  subjection  of  a (111)  grown  sample  to  mechanical 
tension  showed  a shift  in  the  excitonic  absorption  peak  of 
the  MQW  to  lower  energies,  depending  on  how  much  pressure  was 
introduced,  we  showed  that  an  application  of  a low  bias 
would  produce  a shift  of  24  mev  at  300  K,  and  53  meV  at  13  K. 
This  discrepancy  in  the  value  of  shift  between  the  two 
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temperatures  is  probably  due  to  the  added  tension  because  of 
the  shrinkage  of  the  adhesive  used  to  place  the  sanple  on  the 
test  setup.  Since  the  application  of  strain  can  produce  a 
shift  in  the  cibsorption  peak  , this  feature  can  be  used  to 
operate  this  sanple  as  an  electrically  controlled  switch. 

Fincd.ly,  the  value  of  the  nonlinear  refractive  index  was 
measured,  with  the  help  of  a hybrid  Mach-Zender 
interferometer  setup,  using  CW  TE  lasers  operating  at  980  and 
992.2  nm.  The  980  nm  operating  laser  is  112  meV  and  140  meV 
below  the  Cl-HHl  transitions  of  samples  SI  and  S2, 
respectively,  vrtiile  for  the  992.2  nm  laser  this  separation  is 
127  meV  and  167  meV,  respectively. 

The  relative  phase  shift  of  the  interference  pattern  due 
to  the  increase  of  optical  radiation  was  recorded,  and  the 
loss  coefficient  was  measured  by  the  cut-back  method.  These 
measurements  yielded  a higher  value  of  n2  for  (111)  sanples 
ccxnpared  with  the  (100)  sanples.  The  figure  of  merit  was  70% 
higher  for  SI (111)  than  SI (100),  such  improvement  was  not 
present  for  the  S2  samples,  probably  because  the  screening 
effect  does  not  play  an  ixportant  role  in  the  latter  case. 

Because  of  the  lack  of  equipment,  it  was  not  feasible  to 
perform  this  same  experiment  with  pulsed  laser  sources  in 
order  to  determine  only  the  electronic  part  of  the  nonlinear 
refractive  index.  Eliminating  thermal  effects,  we  expect 
this  value  to  be  higher  also  in  the  case  of  (111)  grown 
structures.  This  can  be  explained  by  the  fact  that  there 
are  two  factors  contributing  to  the  increase  in  the  value  of 


n2eiectronic(the  sign  of  vrtiich  is  negative).  These  factors  are 
the  electronic  effect  working  to  increase  the  bandgap 
reducing  the  refractive  index,  and  the  screening  of  the 
strain-induced  electric  field  which  has  the  same  effect  on 
the  refractive  index.  While  the  first  factor  is  present  in 
(100)  grown  structures,  the  second  effect  is  completely 
absent. 

This  research  showed  that  the  (111)  grown  with  a high 
In  percentage  (10%)  should  be  a better  choice  than  their 
(100)  counterparts  vdien  used  in  the  design  and  operation  of 
nonlinear-based  optical  devices. 


APPENDIX  A 

CALCULATICXI  OF  IHE  Cl-HHl&LHl 
TRANSITION  ENERGIES  IN  (100)  SAMPLES 


The  calculation  of  the  value  of  the  C1-HH1,LH1 
transition  in  the  InGaAs/AlGaAs  MQWs  is  performed  using  the 
Kronig-Penney  model.  The  effect  of  strain  on  these 
transitions  is  also  included. 

For  the  case  of  an  infinite  periodic  one-dimensional 
square-vrell  potential  (Fig.  A.l)  it  is  possible  to  arrive  at 
an  exact  solution  of  the  Schr  odinger  equation.  Kronig  and 
Penney  investigated  the  solution,  and  although  it  is  related 
to  a somewhat  idealized  periodic  potential  which  is  only  a 
crude  approximation  to  that  found  in  the  actual  crystal,  it 
is  nevertheless  very  useful  because  it  serves  to  illustrate 
in  a most  explicit  way  many  of  the  inportant  characteristic 
features  of  the  quantum  behavior  of  electrons  in  periodic 
lattices.  The  wave  equation  associated  with  this  model  may 
be  calculated  in  the  one-electron  approximation  by  solving 
Schr  odinger 's  equation 


+ ^ (E-V(x)tp(x)  ) = 0 


h 


(A.l) 
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for  a single  electron  in  the  periodic  potential  V(x).  Since 
the  wave  functions  must  have  the  Bloch  form,  it  is  expected 
that 

Tj)(x)  = e^  u(x)  (A.2) 


Substituting  Eg.  A.2  in  Eg.  A.l,  it  is  found  that  the 
function  u(x)  must  satisfy 


+ 2iJc  ^ - 
dx 


^ 2my(x) 


u(x)  = 0 


(A.  3) 


V(x) 

A 


0 

Fig.  A.l  Ideal  periodic  sguare  well  potential. 


v^iere 


a = 


(A.4) 


For  the  potential  in  Fig.  A.l,  it  is  found  that 


103 


(A.5) 


, . du , ( 2 ^2\ 

Y + 2xk  - (k  - p ) U2(x)  = 0 (-b  < X < 0) 


(A.  6) 


where  ui(x)  represents  the  value  of  u(x)  in  the  interval 
(0<x<a)  and  U2(x)  represnts  the  value  of  u(x)  in  (-b<x<0), 
and  vAiere  p is  an  imaginary  value  expressed  by 


where  A,  B,  C,  and  D are  arbitrary  constants. 

The  requirement  of  continuity  for  the  wave  function  and  its 
derivative  at  x = a and  x = -b  demands  that  the  functions 
u(x)  satisfy  the  same  conditions.  Using  these  boundary 
conditions,  the  resulting  equation  is 


(A.  7) 


The  solution  for  these  differential  equations  are 


Ui(x)  = A ei(  a-k)x  + B Q-i{  a+k)x  (0<X<a) 


(A.  8) 


U2(x)  = C ei(  + D e-i(  P+^)x  (-b<x<0) 


(A.9) 


2oc/ 


sinh(Yb)  sin(oa)  + cosh(vb)  cos(oa)  = cos  k(a+b) 


(A.  10) 
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vAiere 

P = iy  (A.ll) 


Solving  Eq.A.lO  will  yield  the  energy  values  of  the  needed 
transitions . 

Furthermore,  the  application  of  biaxial  stress  in  the 
plane  of  growth  can  be  pictured  as  a hydrostatic  stress  and  a 
uniaxial  stress  along  the  direction  perpendicular  to  the 
growth  plane.  The  hydrostatic  ccxiponent  shifts  the  center  of 
gravity  of  the  J = 3/2  multiplet  and  the  J = 1/2  multiplet 
relative  to  the  conduction  band.  The  uniaxial  strain 
component  splits  the  normally  degenerate  J = 3/2,  mj  = 1/2 
and  J = 3/2,  mj  = 3/2  valence  band  edges.  For  layers  under 
compression,  mj  = 3/2  holes  become  heavy  in  the  growth 
direction,  and  light  in  the  layer  plane,  and  vice-versa  for 
the  mj  = 1/2  holes. 

with  the  superposition  of  the  hydrostatic  and  uniaxial 
stress  components,  shifts  in  the  bandgap  can  be  writen  as 


AEo(1) 

AEo(2) 


-2  a 


-2  a 


Cii+  2 

. C,,  /J 

Cjit  2 C 

Cii  /J 


£ 


E 


(A.  12) 
(A.  13) 
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Where  AEq(1)  , AEq(2)  ar  the  shifts  of  the  heavy  hole,  and  the 

light  hole  valence  bands,  respectively.  e is  the  biaxial 
strain  parallel  to  [100]  and  [010]  directions,  taken  here  to 
be  the  mismatch  Aa,  a and  b are  the  hydrostatic  and  the  shear 
deformation  potentials,  respectively. 

The  vcdues  of  a,  b,  Cll  and  C12  for  different  InxGai_xAs 
conpositions  are  calculated  by  interpolating  the  binary 
values  for  GaAs  and  InAs,  using  Vegard's  law. 


Program  to  Calculate  the  C1-HH1,LH1 
Transition  Energy 
In  (100)  Grown  San{)les 

*****  Deriving  the  Parameters  of  AlGaAs  an  InGaAs 

clc 

clear 

M=9.1095e-31; 

x=0.1; 

y=0.2; 

i=l; 

j=i; 

lc=l; 

hbl  = 6.58219e-16; 
hb2  = 1.05459e-34; 

aGaAs  = 5.6533; 
alnAs  = 6.0584; 
aAlAs  = 5.6611; 

aInGaAs  = x*aInAs  +(l-x)*aGaAs; 
aAlGaAs  = y*aAlAs  +(l-y)*aGaAs; 


CllGaAs  = 11.88; 
C12GaAs  = 5.38; 


CllInAs  = 8.329; 
C12lnAs  = 4.526; 


***** 


CllInGaAs  = x*CllInAs  + (l-x)*CllGaAs; 
C12lnGaAs  = x*C12lnAs  + (l-x)*C12GaAs; 

AGaAs  = -9.8; 

AlnAs  = -5; 

BGaAs  = -1.76; 

BlnAs  = -1.8; 


AInGaAs  = x*AInAs  +(l-x)*AGaAs; 

BInGaAs  = x*BInAs  +(l-x)*BGaAs; 

EglnGaAs  = 0.36  + 0.505  *(l-x)  + 0.555*( (l-x)"2) 
EgAlGaAs  = 1.424  + 1.266  * y + 0.26  *(y"2); 

mcInGaAs  = M*(x*0.023  + (l-x)*0.067) ; 
itcAlGaAs  = M*(0.0665  + 0.0835  * y); 

mhhlnGaAs  = M*(x*0.34+(l-x)*0.92); 
mhhAlGaAs  = M*(y*0.4+(l-y)*0.35) ; 

mlhlnGaAs  = M*(x*0.027  + (l-x)*0.08); 
mlhAlGaAs  = M*(0.18*y  + (l-y)*0.08); 


*****  The  Energy  Shift  Due  to  Strain  ***** 


Aa  = (aInGaAs-aAlGaAs)/aAlGaAs; 

cJEhhl  = -2*AInGaAs*((ClllnGaAs-C12lnGaAs) /CllInGaAs) 


dEhh2  = BInGaAs*((CllInGaAs+2*C12lnGaAs)/CllInGaAs); 
dEhh  = (dEhhl  + dEhh2)*  Aa; 

dElhl  = -2*AInGaAs*((CllInGaAs-C12InGaAs)/CllInGaAs) 
dElh2  = -BInGaAs*((CllInGaAs+2*C12lnGaAs)/CllInGaAs) 
dElh  = (dElhl  + dElh2)*  Aa; 

*****  Kronig-Penney  Model  ***** 

lb=150e-10; 

lz=100e-10; 

Ee  = hbl*hb2*((pi/(l2+27.5e-10))"2)/(2*ncInGaAs) 

Evh  = hbl*hb2*((pi/(lz+1.5E-10))"2)/(2*mhhInGaAs) 

dEv  = 0.4  * (1.247*y  + 1.5*x  - 0.4*(x"2)); 
dEc  = 0.6  * (1.247*y  + 1.5*x  - 0.4*(x"2)); 

Eg  = zeros (1, 2000 ) ; 

Evhh  = zeros ( 1,2000) ; 

Evlh  - zeros (1, 2000 ) ; 

for  E=le-3:le-3:dEv 

alphal=sqrt ( 2 *E*ncInGaAs / ( hbl *hb2 ) ) ; 

gaiimal=sqrt  ( 2 * ( dEc-E ) *incAlGaAs  / ( hbl  *hb2 ) ) ; 

rl=  ( ( gamroal ) ^2-  ( alphal ) "2 ) / ( 2*alphal  *gainnal ) ; 


r2  = (sinh(ganinal*lt))*sin(alphal*lz) ); 
r3  = (cosh(ganinal*lb)*(x>s(alphal*lz) ); 
r=rl*r2+r3; 
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if  sign(r)<=0 

Ec(l,i)=E; 

i=i+l; 

end 

alpha2=sqrt  ( 2 *E*mhhInGaAs  / ( hbl  *hb2 ) ) ; 
gainma2=sqrt  ( 2 * ( dEv-E ) *mhhAlGaAs  / ( hbl  *hb2 ) ) ; 


tl=  ( ( gamma2 ) ^2- ( alpha2  )^2)/(  2*alpha2*ganma2 ) ; 
t2  = (sinh(ganina2*lb)*sin(alpha2*lz) ); 
t3  = (cosh(ganma2*lb)*cx>s(alpha2*lz) ); 
t=tl*t2+t3; 


if  sign(t)<=0 
Evhh(l,j)=E; 

j=j+i; 

end 

alpha3=sqrt(2*E*inlhInGaAs/(hbl*hb2) ) ; 
ganma3=sqrt  ( 2 * ( dEv-E ) *mlhAlGaAs  / ( hbl  *hb2 ) ) ; 


sl=  ( (gainna3 ) "2-  ( alpha3 ) "2 ) / ( 2*alpha3*gamna3 ) ; 


s2  = (sinh(ganiiia3*lb)*sin(alpha3*lz) ) ; 
s3  = (cosh(gainna3*lb)*cos(alpha3*lz) ); 
s=sl*s2+s3; 
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if  sign(s)<0 
Evlh(l,k)=E; 
k=k+l; 
end 

end 


*****  Cl-HHI  and  Cl-LHl  Transition  Energies 

Ellh  = EglnGaAs  + Ec(l,l)  + Evhh(l,l); 

Elll  = EglnGaAs  + Ec(l,l)  + Evlh(l,l); 

Ellh  = Ellh  + dEhh 
Elll  = Elll  + dElh 
Ellhw  = 1.24/Ellh 
Elllw  = 1.24/Elll 
Ec(l,l) 

Evhh(l,l) 


***** 


APPENDIX  B 

CALCULATICXJ  OF  IHE  Cl-HHl&LHl 
TRANSITION  ENERGIES  IN  (111)  SAMPLES 

The  calculation  of  the  value  of  the  C1-HH1,LH1 
transition  in  the  InGaAs/AlGaAs  MQWs  is  done  using  the 
Kronig-Penney  model.  The  effect  that  strain  has  on  these 
transitions  is  also  included.  In  the  case  of  (111)  grown 
layers,  the  energy  shift  due  to  strain  is  larger  than  that 
in  the  (100)  counterparts.  Furthermore,  the  presence  of  a 
strain- induced  electric  field  in  the  strained  layers  would 
reduce  the  transitions  energies.  These  different  factors  cure 
incorporated  in  the  following  program  to  calculate  the  values 
of  the  needed  transitions. 
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Program  to  Calculate  the  Cl-HH1,LH1 
Transition  Energy 
In  (111)  Grovm  Sanples 
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*****  Deriving  the  Parameters  of  AlGaAs  an  InGaAs 


clc 

clear 

M=9.1095e-31; 

x=0.1; 

y=0.15; 

i=i; 

j=i; 

k=l; 

hbl  = 6.58219e-16; 
hb2  = 1.05459e-34; 
aGaAs  — 5 • 6533; 
alnAs  = 6.0584; 
aAlAs  = 5.6611; 

aInGaAs  = x*aInAs  +(l-x)*aGaAs; 
aAlGaAs  = y*aAlAs  +(l-y)*aGaAs; 
o = (alnGaAs-aAlGaAs)/aAlGaAs; 
fa=0; 
fb=0; 


***** 


eGaAs  = 13.1; 
eAlAs  = 10.1; 


elnAs  = 15.15; 


eAlGaAs 

= y*eAlAs  + (l-y)*eGaAs; 

eInGaAs 

= X * elnAs  + (1  - X)  * eGaAs; 

el4GaAs 

VO 

t 

1 

II 

el4lnAs 

= -.045; 

el4lnGaAs  = x * el4InAs  + (1  - x)  * el4GaAs; 


CllGaAs 

= 11.88; 

C12GaAs 

= 5.38; 

CllInAs 

= 8.329; 

C12lnAs 

= 4.526; 

C44GaAs 

= 5.92; 

C44lnAs 

= 3.96; 

CllInGaAs  = x*CllInAs  + (l-x)*CllGaAs; 
C12lnGaAs  = x*C12lnAs  + (l-x)*C12GaAs; 
C44lnGaAs  = x * C44lnAs  + (1  - x)  * C44GaAs; 


AGaAs  = ■ 

1 

VO 

• 

00 

AlnAs  = ■ 

-5* 

BGaAs  = ■ 

-1.76; 

BlnAs  = • 

-1.8; 

AInGaAs  = x*AInAs  +(l-x)*AGaAs; 
BInGaAs  = x*BInAs  +(l-x)*BGaAs; 


acGaAs  = -18.3; 
acInAs  = -11.7; 
avGaAs  = -4.8; 
avInAs  = -5.2; 
dGaAs  = -4.4; 
dInAs  = -4.6; 

acInGaAs  = x*acInAs  +(l-x)*acGaAs; 
avInGaAs  = x*avInAs  +(l-x)*avGaAs; 
dInGaAs  = x*dInAs  +(l-x)*dGaAs; 

mcInGaAs  = M*(x*0.023  + (l-x)*0.067) ; 
mcAlGaAs  = M*(0.0665  + 0.083  * y); 

mhhlnGaAs  = M*0.92; 

mhhAlGaAs  = M*(l.ll*y+0.92*(l-y) ); 

ndhlnGaAs  = M*(x*0.026  + (l-x)*0.07); 
mlhAlGaAs  = M*(0.14*y  + (l-y)*0.07); 

EglnGaAs  = 0.36  + 0.505  *(l-x)  + 0.555*( (l-x)''2) ; 
EgAlGaAs  = 1.424  + 1.266  * y + 0.26  *(y"2); 

Ee  = hbl*hb2*((pi/(lz+27.5e-10))"2)/(2*ncInGaAs); 
Evh  = hbl*hb2*((pi/(lz+1.5E-10))"2)/(2*mhhInGaAs) 


dEv  = 0.4  * (1.247*y  + 1.5*x  - 0.4*(x"2)); 
dEc  = 0.6  * (1.247*y  + 1.5*x  - 0.4*(x"2)); 

Ec  = zeros (1,2000); 

E\^  = zeros  (1, 2000 ) ; 

Evlh  = zeros (1, 2000 ) ; 


*****  Amount  of  Electric  Field  Inside  the  MQW  Region  ***** 

Aa  = (CllInGaAs+  2 * C12lnGaAs)  / (4  * C44lnGaAs); 

Ab  = (CllGaAs  + 2 * C12GaAs)  / (4  * C44GaAs); 

B = (CllInGaAs  + 2 * C12lnGaAs)  / (CllGaAs  + 2 * C12GaAs); 

exxa  = ((aGaAs  / (aInGaAs))  - 1)  / ((1+Aa)  + ((1+Ab)  *B 

* (Iz  /lb))); 

exxb  = -(Iz  / lb)  * B * exxa; 

exya  = -2  * Aa  * exxa; 

exyb  = -2  * Ab  * exxb; 

fa=  2 * exya  * el4InGaAs  / (8.8542E-12  * eInGaAs) 
fb=  2 * exyb  * el4GaAs  / (8.8542E-12  * eAlGaAs); 


*****  The  Energy  Shift  Due  to  Strain  ***** 
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DeltaEc  = 3*acInGaAs*exxa; 

DeltaEhh  = 3*avInGaAs*exxa  + 1.73205*dInGaAs*exya; 
DeltaElh  = 3*avInGaAs*exxa  - 1.73205*dInGaAs*exya; 


*****  Kronig-Penney  Model  ***** 

lb=178e-10; 

lz=80e-10 

for  E=le-3:le-3:dEv 
alphal=sqrt(2*E*mcInGaAs/(hbl*hb2 ) ) ; 
ganinal=sgrt  ( 2 * ( dEc-E ) *mcAlGaAs  / ( hbl  *hb2 ) ) ; 

rl=  ( (gammal)''2-(alphal)^2)/(2*alphal*gainmal) ; 
r2  = (sinh(ganmal*lb)*sin(alphal*lz) ); 
r3  = (cosh(gaiiinal*lb)*cos(alphal*lz) ); 
r=rl*r2+r3; 

if  sign(r)<=0 
Ec(l,i)=E; 
i=i+l ; 


end 


cLLpha2=sqrt  ( 2 *E*mhhInGaAs  / ( hbl  *hb2 ) ) ; 
ganina2=sqrt  ( 2 * ( dEv-E ) *inhhAlGaAs  / ( hbl  *hb2 ) ) ; 
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tl=  ((ganina2)''2-(alpha2)''2)/(2*alpha2*gaiiina2); 
t2  = (sinh(gamma2*lb)*sin(alpha2*lz) ); 
t3  = (cosh(gainna2*lb)*cos(alpha2*lz) ) ; 
t=tl*t2+t3; 

if  sign(t)<=0 
Evhh(l, j)=E; 

j=j+i; 

end 

alpha3=sgrt ( 2 *E*mlhInGaAs / ( hbl *hb2 ) ) ; 
gamna3=sqrt  ( 2 * ( dEv-E ) *mlhAlGaAs  / ( hbl  *hb2 ) ) ; 

sl=  ( (ganma3)"2-(alpha3)"2)/(2*alpha3*gainma3); 
s2  = (sinh(ganma3*lb)*sin(alpha3*lz) ); 

S3  = (cosh(ganina3*lb)*cos(alpha3*lz)); 
s=sl*s2+s3; 

if  sign(s)<0 
Evlh(l,k)=E; 
k=k+l; 
end 


end 
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*****  Electric  Field  Effect  on  C1-HH1,LH1  Transitions  ***** 


Lzc  = 6.626e-34*.5/(sqrt(2*mcInGaAs*Ec(l,l)*1.6e-19) ) ; 
Lzhh  = 6.626e-34*.5/(sqrt(2*mhhInGaAs*Evhh(l,l)*1.6e-19) ) ; 
Lzlh  = 6.626e-34*.5/(sqrt(2*mlhInGaAs*Evlh(l,l)*1.6e-19) ); 

Flc  = Ec(l,l)/Lzc; 

Flhh  = Evhh( 1,1) /Lzhh; 

Fllh  = Evlh( 1,1) /Lzlh; 

fc  = abs(fa/Flc); 
fvdih  = abs(faZFlhh) ; 
fvlh  = abs(fa/Fllh); 

WC=Ec(l,l)*(0. 99969+0. 0144*fc-0. 014364* (fc"2)+0.00042568 
*(fc"3)- 

5.58e-6*(fc''4) ); 

Whh  = Evhh(l,l)*(0. 99969+0. 0144*fvhh-0.014364*(fvhh"2)  + 

0 . 00042568* ( fvhh"3 ) -5 . 58e-6* ( fvhh"4 )) ; 

Wlh  = Evlh(l,l)  * (0. 99969+0. 0144*fvlh-0.014364*(fvlh^2)+ 

0 . 00042568*  ( fvlh''  3 ) -5 . 58e-6*  ( fvlh"4 ) ) ; 


*****  Cl-HHl  and  Cl-LHl  Transition  Energies 


***** 


Ellh  = EglnGaAs  + Wc  + Whh  + DeltaEc  + DeltaEhh 
Elll  = EglnGaAs  + Wc  + Wlh  + DeltaEc  + DeltaElh 


Ellhw  = 1.24/Ellh; 
Elllw  = 1.24/Elll; 


APPENDIX  C 

CALCULATION  OF  THE  REFRACTIVE  INDICES 
FOR  IHE  LAXERS  IN  STRUCTURES  SI  AND  S2 


Program  to  Calculate  the  Refractive  Indices 
Of  the  Different  Layers  in  Structures  SI  & S2 


clc 

clear 

clg 

itFzeros  (4,1); 
nGaAs  = zeros(l,500) ; 
nInGaAs  = zeros(l,500) ; 
Wavelength  = zeros ( 1 , 500 ) ; 
Energyl  = zeros (1, 500 ) ; 
Wavelengthl  = zeros (1,500); 


y=0.2; 

EgAlGaAs=l . 424+1 . 266*y+0 . 26*y^'2 ; 
EoAlGaAs=3 . 65+0 . 871*y+0 . 179*y"2 ; 
EdAlGaAs=36 . 1-2 .45*y ; 
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x=0.1; 

EgInGaAs=0 . 36+0 . 505* ( 1-x ) +0 . 555* ( ( 1-x ) " 2 ) ; 

EoInGaAs=3 .65-2 . 15*x; 

EdInGaAs=36 . 1-19 . 9*x; 

EgGaAs=1.42; 

EoGaAs=3.65; 

EdGaAs=36.1; 

il=(EgAlGaAs+0.3)*500  + 1; 
nAlGaAs  = zeros(l,il); 

Energyl  = zeros(l,il); 

Wavelengthl=zeros ( 1 , il ) ; 
t=l; 

for  E=0.2:0.002:EgAlGaAs+0.5 
Energy l(l,t)=E; 

Wavelengthl ( 1 , t ) =1 . 24 /E ; 
p=l+( (E/EoAlGaAs ) "2 ) ; 

q=(E^4/(2*(EoAlGaAs"2)*( (EQAlGaAs"2)-(EgAlGaAs"2) ) ) ) ; 
v=log((2*(EoAlGaAs''2)-(EgAlGaAs''2)-(E"2) )/( (EgAlGaAs"2)- 
E-2)); 

nAlGaAs  ( 1 , t ) = ( 1+ ( EcJAlGaAs /EoAlGaAs ) * ( p+q*  V ) ) " 0 . 5 ; 

t=t+l; 

end 


i2=(EgInGaAs+0.3)*500  + 1; 


nInGaAs  = zeros(l,i2); 

Energy2  = zeros(l,i2); 

Wavelength2=zeros ( 1 , i2 ) ; 
t=l; 

for  E=0.2:0.002:EgInGaAs+0.5 
Energy2(l,t)=E; 

Wavelength2 ( 1 , t ) =1 . 24  /E ; 
p=l+((E/EoInGaAs)"2); 

q=(E^4/(2*(EoInGaAs"2)*(  (EoInGaAs''2)-(EgInGaAs^2) ) ) ) ; 
v=log( (2*(EoInGaAs"2)-(EgInGaAs"2)-(E"2) )/( (EgInGaAs"2)- 
E-2)); 

nInGaAs  ( 1 , t ) = ( 1+ ( EdInGaAs /EoInGaAs ) * ( p+q*v ) ) '^  0 . 5 ; 

t=t+l; 

end 

i3=(EgGaAs+0.3)*500  + 1; 
nGaAs  = zeros(l,i3); 

Energy3  = zeros (1,13); 

Wavelength3=zeros (1,13); 
t=l; 

for  E=0.2:0.002:EgGaAs+1.5 
p=l+( (E/EoGaAs)"2); 

q=(E^4/(2*(EoGaAs"2)*( (EoGaAs^2)-(EgGaAs"2) ) ) ); 
v=log( (2*(EoGaAs^2)-(EgGaAs"2)-(E^2) )/( (EgGaAs"2)-E^2) ) ; 
nGaAs  ( 1 , t ) = ( 1+ ( EdGaAs /EoGaAs ) * ( p+q* V ) ) '^  0 . 5 ; 

Energy 3 (l,t)=E; 

Wavelength3 ( 1 , t ) =1 . 24 /E ; 
t=t+l; 


end 
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*****  Plots  of  Refractive  Indices  of  ***** 

*****  GaAs  AlGaAs  and  InGaAs  vs.  Energy  ***** 

m(l,l)=0.2;m(2,l)=1.6;m(3,l)=3.2;m(4,l)=3.8; 

axis(m) 

plot  (Energy  1 ( 1 , 1 : il ),  nAlGaAs  ( 1 , 1 ; il ),  Energy2  ( 1 , 1 : i2 ),  nInGaAs  ( 
l,l:i2) , ' — ' , Energy 3, nGaAs, ) 
title ( 'Refractive  Index  vs  Energy') 
xlabel( ' Energy (eV) ' ) 

ylabel(  'n(InGaAs-top,GaAs-middle,AlGaAs-bottom) ' ) 
grid 

text (1,3. 7, 'InGaAs  — ') 
text(l,3.65, 'GaAs  -. ' ) 
text( 1,3.6, 'AlGaAs  -') 


APPENDIX  D 

CONDITIONS  FOR  SINGLE  MODE  OPERATION 


The  following  program  uses  the  multilayer  matrix  method 
to  find  the  propagation  constants  of  the  fundamental  mode  for 
the  san^les  SI  cind  S2,  making  sure  that  the  choice  of  layer 
thicknesses  and  conpositions  is  suitable  for  single  lateral 
mode  operation.  Furthermore,  the  mode  field  of  these 
structures  is  plotted,  showing  how  the  optical  field  is 
distributed  inside  each  structure.  The  refractive  indices 
for  the  different  layers  are  derived  in  Appendix  C. 

The  last  part  of  the  program  deals  with  the  design  of  a 
single  mode  ridge  waveguide,  using  the  effective  index 
method.  If  a certain  waveguide  width  is  needed,  we  consider 
the  structure  to  be  etched  to  a certain  depth.  Then  we  find 
the  effective  index  in  such  structure.  Using  some  iteration, 
the  amount  of  etch  (or  the  width  of  the  waveguide)  can  be 
determined  to  insure  single  mode  operation. 
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Program  to  Obtain  the  Propagation  Constant  & 
Mode  Field  in  Proposed  Structures 
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clc 

clg 

clear 


*****  Values  of  the  Refractive  Indices  & ***** 
*****  Thicknesses  for  Different  Layers  ***** 


nl=l; 

n2=3.53;h2=0.05; 

n3=3.419;h3=l; 

n4=3.438;h4=0.25;h41=0.2; 

n5=3.4578;h5=0.25; 

n6=3.438;h6=.25; 

n7=3.419; 

X=0.9538; 

*****  Constants  and  Intialiazation  of  Variables 

j=sqrt(-l); 

A=0; 

B=0; 

C=0; 

E>=0; 

E=0; 

P=0; 

G=0; 


***** 
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H=0; 

1=0; 

J=0; 

el=0; 

W=0; 

k0=2*pi/  X,; 
cnt2=0;z=0; 

ul=(n5-n3-0.0002)*10000; 
u=zeros(ul,l) ; 
v=zeros (4,1); 
x=zeros(5,l) ; 
y=zeros(5,l) ; 
cnt3=0; 
cnt4=l ; 


ol=0 

o2=0 

o3=0 


o4=0 


r=zeros (1,1); 


s=0; 


fl=zeros(ul,l); 

f2=zeros(ul,l); 


c=zeros(2,2) ; 
iiF  zeros (2,2); 
m2=zeros(2,2) ; 
m3=zeros (2,2); 
m4=zeros(2,2); 


in41=zeros(2,2); 
m5=zeros(2,2); 
m6=zeros(2,2) ; 
t=zeros(2,2) ; 
cnt=2 ; 

*****  Multilayer  Matrix  Method  ***** 

for  I^n5-0.0001:-0.0001:n3+0.0001 

u(cnt4,l)=N; 

i2=n2"2-N''2; 

i3=n3"2-N''2; 

i4=n4"2-N"2; 

i5=n5"2-N"2; 

i6=n6"2-N"2; 

ka2=k0*sqrt(abs(i2) ) ; 

if  i2<0 

ka2=ka2*j; 

end 

ka3=k0*sqrt(abs(i3) ) ; 

if  i3<0 

ka3=j*ka3; 

end 


ka4=k0*sqrt ( abs ( i4 ) ) ; 
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if  i4<0 

ka4=ka4*j; 

end 


kci5=k0*sqrt  ( abs  ( i5 ) ) ; 

if  i5<0 

ka5=ka5*j; 

end 

ka6=k0*sqrt ( abs ( i6 ) ) ; 

if  i6<0 

ka6=ka6*j; 

end 

m2(l,l)=cos(ka2*h2) ; 
m3  ( 1 , 1 )=c»s ( ka3*h3 ) ; 
m4(l,l)=<x>s(ka4*h4) ; 
m41 (1,1 )=cos (ka4*h41 ) ; 
m5  ( 1 , 1 )=c»s ( ka5*h5 ) ; 
m6  ( 1 , 1 ) =cx5S  ( ka6*h6 ) ; 

m2(2,2)=m2(l,l); 

m3(2,2)=m3(l,l); 

m4(2,2)=m4(l,l); 
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m41(2,2)=m41(l,l); 

m5(2,2)=m5(l,l); 

m6(2,2)=m6(l,l) ; 

m2(l,2)=( j*sin(ka2*h2) )/ka2; 
m2 ( 2 , 1 ) = j *ka2 *sin ( ka2 *h2 ) ; 

m3 ( 1 , 2 )= ( j*sin(ka3*h3 ) ) /ka3 ; 
m3(2,l)=j*ka3*sin(ka3*h3) ; 

m4(l,2)=( j*sin(ka4*h4) )/ka4; 
m4 { 2 , 1 ) = j *ka4 *sin ( ka4 *h4 ) ; 

m41(l,2)=( j*sin(ka4*h41) )/ka4; 
m41(2,l)=j*ka4*sin(ka4*h41) ; 

m5(l,2)=( j*sin(ka5*h5) )/ka5; 
m5 ( 2 , 1 )= j *ka5*sin (ka5*h5 ) ; 

m6(l,2)=( j*sin(ka6*h6) )/ka6; 
m6(2,l)=j*ka6*sin(ka6*h6) ; 

ro=m6*m5*m4*m3*m2 ; 
t=m6*m5*m41; 

m(l,2)=-m(l,2)*j; 

m(2,l)=-m(2,l)*j; 
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t(l,2)=-t(l,2)*j? 

t(2,i)=-t(2,i)*j; 

gc=kO*sqrt(N''2-nl''2) ; 
gs=kO*sqirt(N^2-n7^2) ; 

fl(cnt-l,l)=atan(  (gs*m(l,l)+gc*m(2,2)  )-m(2,l)+gc*gs*m(l,2) ) ; 
f 2 (cnt-1 , 1 )=atan( (gs*t( 1 , 1 )+gc*t(2 ,2 ) )-t (2 , 1 )+gc*gs*t ( 1, 2 ) ) ; 

if  cnt>2 

ol=sign( f 1 ( cnt-2 , 1 ) ) ; 
o2=sign(fl(cnt-l,l) ) ; 
if  (ol-K>2)=0 
cnt2=cnt2+l ; 
x(cnt2)=N; 
beta2=kO*N; 
end 

o3=sign( f2 (cnt-2 , 1 ) ) ; 
o4=sign(f2(cnt-l,l) ) ; 

if  (o34o4)=0 
cnt3=cnt3+l ; 
y(cnt3)=N; 
betal=kO*N; 
end 

cnt4=cnt4+l ; 
end 


cnt=cnt+l ; 
end 
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*****  Using  the  Effective  Index  Method  to  ***** 
*****  Calculate  the  Width  of  the  Ridge  Guide  ***** 


cnt6=l ; 

for  d=0.1:0.1:8 

r(l,l)=k0*d*sqrt(x(l,l)"2-y(l,l)"2); 

if  r(l,l)>=pi  & cnt6=l 

s=d; 

MaxWidth jmRB 
cnt6=cnt6+l ; 

end 

end 


V^round(  MaxWidth jm)j 
if  W-MaxWidth jm  >0 
W=W-1; 
end 

if  abs  (W-Ma:rt^idth _^im)<le-4 

W=W-0.1; 

end 

W 


***  Using  tte  Multilayer  Matrix  Method  to  Determine 
***  The  Mode  Field  and  Effective  Index  *** 

***  of  Etched  Structure  *** 


*** 
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cnt2=0; 

u2=(x(l,l)-y(l,l)-0.0001)*100000; 

f3=zeros(u2,l); 

g=zeros(u2,l); 

cnt=l; 

for  i=y(l,l):0.00001:x(l,l)-0. 00001 

cnt=cnt+l ; 

g(cnt-l,l)=i; 

il=x(l,l)"2-i"2; 

kl=k0*sqrt ( abs ( il ) ) ; 

c(l,l)=cos(kl*(W) ); 

c(l,2)=(j*sin(kl*(W)))/kl; 

c(2,2)=c(l,l); 

c(2,l)=j*kl*sin(kl*(W)); 
c(l,2)=-c(l,2)*j; 
c(2,l)=^(2,l)*j; 
gc=k0*sqrt(i^2-y(l,l)^2) ; 
gs=k0*sqrt(i^2-y(l,l)^2) ; 

f3(cnt-l,l)=atan((gs*c(l,l)+gc*c(2,2))-c(2,l)+gc*gs*c(l,2)); 
if  cnt>2 

ol=sign( f 3 (cnt-2 , 1 ) ) ; 
o2=sign(f3(cnt-l,l) ) ; 

if  (ol+o2)=0 
z=i; 


end 
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end 

end 

El=zeros ( ( h2+h3+h4+h5+h6+4 ) *100+1 , 1 ) ; 
E2=zeros ( (ti41+h5+h6+4 ) *100+1 , 1 ) ; 
E3=zeros ( ( (Wf4 ) *100 )+l , 1 ) ; 

dl=zeros ( ( h2+h3+h4+h5+h6+4 ) *100+1 , 1 ) ; 
d2=zeros ( (h41+h5+h6+4 ) *100+1 , 1 ) ; 
d3=zeros ( ( (Wf4 ) *100 )+l , 1 ) ; 

thetal=asin ( X ( 1 , 1 ) /n5 ) ; 
theta2=asin (y ( 1 , 1 ) /n5 ) ; 
theta3=asin(z/x(l,l) ) ; 

Kappal=k0*n5*cos (thetal ) ; 
Kappa2=k0*n5*cos ( theta2 ) ; 

Kappa3=k0  *x ( 1 , 1 ) *c»s ( theta3 ) ; 

Gall=k0*sqrt(x(l,l)^2-nl^2); 
Gal2=k0*sqrt (x( 1 , 1 ) ^2-n2"2 ) ; 

Gal  3=k0*sqrt  ( X ( 1 , 1 ) 2-n3 2 ) ; 
Gal4=k0*sqrt  ( X ( 1 , 1 ) '' 2-n4  " 2 ) ; 
Gal6=k0*sqrt(x(l,l)^2-n6"2); 

Gal  7=k0*sqrt  (X  ( 1 , 1 ) 2-n7 2 ) ; 
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Ga21=k0*sqrt(y(l,l)^2-nl''2) ; 

Ga24=k0*sqrt(y(l,l)^2-n4^2) ; 

Ga26=k0*sqrt(y(l,l)''2-n6^2) ; 

Ga27=k0*sqrt (y ( 1 , 1 ) "2-n7 "2 ) ; 

Ganina3=k0*sqrt(  (x(l,l)*sin(theta3)  )''2-y(l,l)''2) ; 

A=( (0.5*c»s(Kappal*(h5/2) ) )+(KaEpal*0.5*sin(Kappal*(h5/2) )/Ga 
14)); 

B=( (0.5*cos(Kappal*(h5/2)))- 
(Kappal*0 . 5*sin(Kappal* (h5/2 ) ) /Gal4 ) ) ; 

D=0.5*(A*exp(-Gal4*h4)+B*exp(Gal4*h4))+  (0.5*Gal4/Gal3)*(- 
A*exp(-Gal4*h4)+B*exp(Gal4*h4) ) ; 
C=-IH(A*exp(-Gal4*h4)+B*exp(Gal4*h4) ) ; 

F^O . 5* (C*exp( -Gal3*h3 )+D*exp(Gal3*h3 ) )+  ( 0 . 5*Gal3/Gal2 ) * (- 
C*exp ( -Gal 3 *h3 ) +D*exp ( Gal 3 *h3 ) ) ; 
E=-F+(C*exp(-Gal3*h3)+D*exp(Gal3*h3) ) ; 
G=E*exp(-Gal2*h2)+F*exp(Gal2*h2 ) ; 

I=( (0.5*c»s(Kappal*(h5/2) ) )+(Kappal*0.5*sin(Kappal*(h5/2) )/Ga 
16)); 

H=((0.5*cos(Kappal*(h5/2)))- 
( Kappal *0 . 5*sin ( Kappal * ( h5/2 ) ) /Gal 6 ) ) ; 

J=H*exp ( Gal 6 *h6 ) +l *exp ( -Gal 6 *h6 ) ; 


Al=( (0.5*cos(Kappa2*(h5/2) ) )+(Kappa2*0.5*sin(Kappa2*(h5/2) )/G 
a24); 
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Bl=( (0.5*cos(Kappa2*(h5/2) ) )- 
(Kappa2*0.5*sin(Kappa2*(h5/2) )/Ga24) ) 

Il=( (0.5*cos(Kappa2*(h5/2) ) )+(Kappa2*0.5*sin(Kappa2*(h5/2) )/G 
a26)) 

Hl=  ( ( 0 . 5*<X)s  ( Kappa2*  ( h5/2 ) ) ) - 
(Kappa2*0.5*sin(Kappa2*(h5/2) )/Ga26) ) 
Jl=Hl*e:5)(Ga26*h6)+Il*exp(-Ga26*h6) ; 

Gl=Al*exp(-Ga24*h41  )+Bl*exp(Ga24*h41 ) ; 

M=l; 

for  i=-2-h6;0.01:h5+h4+h3+h2+2 
dl(M,l)=i; 


if  i>=0  & i<=h5 

El(M,l)=cos(Kappal*(i-(h5/2) ) ); 
end 

if  i>h5  & i<=  h5+h4 

El(M,l)=A*exp(-Gal4*(i-h5) )+B*exp(Gal4*(i-h5) ); 
end 


if  i>h5+h4  & i<=h3+h4+h5 

El(M,l)=C*exp(-Gal3*(i-h5-h4) )+D*exp(Gal3*(i-h5-h4) ) ; 
end 


if  i>h3+h4+h5  & i<=  h2+h3+h4+h5 
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El(M,l)=E*exp(-Gal2*(i-h5-h4-h3))  + F*exp(Gal2*(i-h5-h4-h3) ) ; 
end 

if  i>h2+h3+h4+h5 

El(M,l)=G*exp(-Gall*(i-h2-h3-h4-h5) ) ; 
end 


if  i<0  & i>=-h6 

El (M, 1 )=I*exp(Gal6*i)+H*exp(-Gal6*i) ; 
end 


if  i<-h6 

El  (M, 1 )=J*exp(Gal7* (i+h6 ) ) ; 

end 

»=M+1; 

end 

M=l; 

for  i=h41+2;-0.01:-h5-h6-2 

d2(M,l)=i; 

if  i<0  & i>=-h6 

E2(M,l)=Il*exp(Ga26*i)+Hl*exp(-Ga26*i); 

end 


if  i<-h6 

E2 (M, 1 )=Jl*exp(Ga27* (i+h6 ) ) ; 
end 
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if  i>=0  & ic=h5 

E2 (M, 1 )=cos (Kappa2* ( i- ( h5/2 ) ) ) ; 
end 


if  i>h5  & i<=  h5+h41 

E2(M,l)=Al*exp(-Ga24*(i-h5) )+Bl*exp(Ga24*(i-h5) ) ; 
end 

if  i>h41+h5 

E2(M,l)=Gl*exp(-Ga21*(i-h41-h5) ) ; 

end 

M=M+1; 

end 

M=l; 

for  i=2+W:-0.01:-2 

d3(M,l)=i; 

if  i>=0  & i<=W 

E3 (M, 1 )=cos ( Kappa3* ( i- (W/2 ) ) ) ; 
end 

if  i>W 

E 3 ( M , 1 ) =cos  ( K^pa3  *W/2 ) *exp  ( -Ganna3  * ( i-W ) ) ; 

end 

if  i<0 

E3 ( M , 1 ) =cos ( Kappa3  *W/ 2 ) *exp ( Ganroa3  *i ) ; 
end 


end 
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*****  Graphs  of  Effective  Indices  of  the  Whole,  **^ 
*****  The  Etched  and  the  Three-Layer  Structures 

v(l,l)=n3;v(2,l)=n5;v(3,l)=-2;v(4,l)=2; 

axis(v) 

subplot (221) 

plot(u,fl(l:ul,l) ) 

title ( 'Effective  Index  for  the  Whole  Structure') 

xlabel ( ' Refractive  Index ' ) 

grid 

subplot (222) 
plot(u,f2(l:ul,l) ) 

title ( 'Effective  Index  for  the  Etched  Structure') 

xlabel ( ' Refractive  Index ' ) 

grid 

v(l,l)=y(l,l);v(2,l)=x(l,l);v(3,l)=-2;v(4,l)=2; 

axis(v) 

suh^lot(223) 

plot(g,f3) 

title ( 'Effective  Index  for  the  Three-layer  Structure') 

xlabel ( ' Refractive  Index ' ) 

grid 


pause 

clg 
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*****  Mode  Fields  in  All  Three  Structures  ***** 


V ( 1 , 1 ) =-h6-2 ; V ( 2 , 1 ) =h5+h4+h3+h2+2 ; v ( 3 , 1 ) =0 ; v( 4 , 1 ) =1 • 5 ; 

axis(v) 

subplot (221) 

plot (dl, El) 

title ( 'Mode  Field  in  Conplete  Structxire' ) 

xlabel(  'Distance  (im) ' ) 

grid 

v(l,l)=-h6-2;v(2,l)=h5+h41+2;v(3,l)=0;v(4,l)=1.5; 

axis(v) 

subplot  (222) 

plot(d2,E2) 

title ( 'Mode  Field  in  Etched  Structure') 

xlabel(  'Distance  (}m) ' ) 

grid 

v(l,l)=-W;v(2,l)=2*W;v(3,l)=0;v(4,l)=1.5; 

axis(v) 

sutplot(223) 

plot(d3,E3) 

title ( 'Mode  Field  in  Ridge  Waveguide' ) 

xlabel( 'Distance  ()an) ' ) 

grid 
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